


Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1991-09 


Evaluation of absolute positioning using the 
Defense Mapping Agency's GASP program 


Bredthauer, Dennis 


Monterey, California. Naval Postgraduate School 
http://hdl.handle.net/10945/28180 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 
| (8 D U DLEY research materiak and institutional publications created by the NPS community. 
il Ser Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


(NA KNOX appointed — and published -- scholarly author. 

>. LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


i A e GANA Te NU CT AS 
TUE CIA AA RE A vue «e AAA k wee e ARS se 
‘ek is wv ahs Ch LAS hit Sh a Rs A ve CNE FE date be? of, iw ne REX 
r CN CITE n 4 7 WE ENT Ya a 
@ as LE wa CE a Ana 
w ' 4 
x 


















































POT a S pare ms RARE 
A HR ss ya o O À y ES 
- ¡EW y + ie is tb ina D y A Ssh 
“a ' ys ry 5 ce eV RE ity NS Leo A a p NA he IRE NC RL AE aa mae ae he > 
o Te an PKR AUX Me ho ah PARENT og HN GA at A aT aay pre ein 
i aA NARAN k ANAL CORA Ad D PA eh SAAR INNS Ab ARO SI PO Kor ee oper, a S 
| | i l AAA Po l CLR fo h T} TST em CT LL 
het À CEE t i E pt E t CON a hes RARO RU Add Y A As deta bys 
Orne + 4° i ES ay, N rae DES AN AE, NES À NET ER à à ENT ve A ar CAGA RE ap > Q = a 
0 0 7 A ARRET A E te 4% A y E OF Atte na, aM TOR CT rec te a Ad i A 
‘ g OO EE ICAO pS dun e D de. AIT A s . ur oy We) te parte oe Ae a er att Pet Pa ty ek, 
g a a4 DA ata vx ea en $ NAS va TEE EMITE ION UD Le a QUES oe e ce Ar a 
E Wat ws SOC ENST ER LA s CAES pa AA ON Hop AU A A plo D ASA e 
3 LL ES D La CE CET res t poyda k Aa. be hla, ) ra as a Aa x eu de 
g i y y H G Û E | Aa 3 UE à me A A Be Lu de mA Si ROUE A RTE EU ak LA de RARE LEL 2] me A 
: A E ea 19 fim Ame ° LE AA A s ` 
i SANS AN CNE SC et SRNR GT RATES EEE AC Sais 
A i S rt ya Oy ? ratte 7 arte 0 Ary a ` Paru À PL Pose: LA LR Lu Li CI D ESC PT Be A PET en 
D QU RE CIN: E AA oil ER RTE a À 
n at À . ead Agaa p n ET LB SOLE TE CT TEE EN STORE LOL TTL LUS rl 
i E re AA de SO ES EA O US Ke D SAN sig R E (Mer DE RR CNE LA O en raies re 
0 SU sSN Ym Ahh ne Oh A wp ave ak UN LN we! GSES ited | AE Le aN, + mr ar es mars Pe 
> a er erry 7 LE ELLE RO PRET dre PSC O ett aCe oR Tree ty AZ WY Wt NS ara Le an eee 
AO H K LATE ui OA is KA AR bé LA Packed he à v BA E COSO RRA RC AS as aR RAT 
$ f k g LE de 0 wt a A A 
a 5 E E AE SA LA a a8 x PAR K y z LACA IA Py wi Dot: rays DONC RE ROUE me ie PLUS A a Se ine pon perte 
y nt os A MI OP RE AAA O A? DU AOS AE RE Arta RCE aed pere gaa . NAAA Rc Phare =a 
a 0 CE A E . LE ` U NES mate A ad LT a « EE b y 
A A A A Des ea A O RR TA TS LRN Stet RNR Te eT acces RER 
A L . A a | r A A AÑ on a GNT y y o = 
DA Enr ER Ar a Hé rea : A d R w ma a A bs Ñ ON Le y ON A POD A ACA AAA 4 pal Ad ba ROSE ent Peer RE TEST A ETS AS re 
g ‘ ' ` g 5 aar ea) Ly <A cs Aa RA VA VE ais w + 4 Md ye ad 3 er 
j ‘ag Sn a A i nr D OCT C E MA aA, A O LE i ue 4 A G bre Py na de va CE 4 », dy yia? RUE en rs De ST CDS Tee CCA 2 BOR TET os robe 
CEA s ' 0 CEE 4% a J $ ' DEA ET sf x LR JR aw A ES y yu de 
ya oy O EC ; RTE ARS 2» D 5 AA APE AE ie i Pe ays Je PNS ENA ee 20m AOS x RIOS ee ane tenn OEE pine rE: Lee RR 
e S A A ~x = à "e 3 CRT RL a 4 a SARA CEA CES QUE ate a 
RL RÉ T A H u a u À t Ge NE b t a RC K : RC ae) a DO oe ae ACER ei Po arty Se ns IO AAN Pi anar Ne he PRES See ENT oe 
5 CR Wie Cees DAC NS A. CET es, AER ROA UTA E 4 reel Ro AO DES let rh page rol did Blast da ma me eur 
A A rei ead ET E 9 i Gg Er a > y A DA 0 DE CITOS pilaa de kkp Ne de ip x ro pc rh arte CT TEE e pr orne per 
1 ' ses 0 a CRE E à CE: te tya AN sy ths 4 RE) DRE CRE % ee B Ld ih GP a, y N deele aa 53d 2 a Pht th ae enpr i per aE o a SS A 
g 0 CIA oak eet ae r IS n a S Y 0 n MER A eS Supe > n Ss CR CES IS CES CAN an ah ees ep dp e CAES Torri A re Seta hl ray qbo” tah A ls 
E ` CRE F ee g T . . e de aa ry 0.4 PENA O ar u row ta Ty y o > A LR H 
tee BARRE NEO LE OCR eet Se 
a “o. g a Y a “me, 5 ni e Sa pm . a : a an y i po j GI 
z g PR ice RE 5 4 a Je + A M a y va A E ne i ro e pais OL A E NE 5 OO CRE ES LE LA LAS ECOLE alee Rie bale bir aster Eee Ve EE dp] i a i bee id rito SCENE or tthe CRT 
A Oem a Ree, he A e A AO Caer EMO AIR COLA LAY AA ee AS Myre eee COTA rs DR Meme 
z Y g g ` . a . ` EA ELE UT A a QC CE rd o oe 
toe A CCS mr ONE A E i Ake | » oe ve a Are ve KEAR Bes Al HE E E MP fe rs + DES SLR CE us PACE a ara de] NO LT" Re. ee no. 
. o, o OR SA A NET O er pa II TN AENA ear IENA Na cal is mes Frc le Le arora Oh: Se eee ARLES dd. lr lei ya a 
: o E O 5 A g : à 5 Y y 0 COR TER A A 
eS Aa ré j ; A He LE Ari res T Pe eae 3 + Pe Pee AS EE sm rindo a td A dy a, bte gr Nabe OPEN a Vte Tr LL et der 
e P À A 3 fe DEL ‘i UI aa n be CAR SO A LE EN A LOUE" CREER TI EN whe kes ETS ac RE Peed ta) le deudas id he e TT 
Pal Si te oe .. q ara ye CA aa LA 0 $47. s A CHALET ONCE RES RAD CET 
S 0 z van n o ` mi} ary PS eon. eae A dy be gee 
0 o o ce Fea uty ot Fo 
L "UE 


RUE TP NN o AMENER o: ELA 
4 ee is ei dea et a ES A TE 
E AE NC CET DORE DEEE br 
A TAO ET | COD A A 


LR pe ad lcd alate Mathie LE LEONE pote ele Cid 
LA Mebane: eel SA Rae oth 
et J ON O Len CARS CONTE EEE ES Heure 
LU 































































































































































































































































La ¡Mas K aan 
DL EP er ALE a ere CRT NE ET rare tr 
tie DCS RER D PRO ENCRES re oer a A TT nt LT AL Enr ar der net de ee et LT 
À CRETE E LEE E y ERE 5 AMA AOS + ae CA a A E E A ee 
ae LE A LS e. CA a Y >» D EL ©] Bo e [TE ae Rv! = a 
o ` O O A TERT) Ke SAN TN FA 4 a A A; ba o As pa ver LEZ veh ariaya à RTS CA ¿e gt AE yes Cd OL EL SE NE Er Pere Fie eee te 
US . O ES AN MOS ALTO TAE A he D ie: UA w LM eu) RE "y Le do CE pére prie PES à gin ‘= 
o L O 1 MEA AAA ET ET casa e y SE re AS IA O A Aer dl anag, Saa goy y Q D LU 478 ln Mg ph Eh] ER a 
G g 0 à A OI AA sis ed VEIS ree RE LCL NS A ÍA 
A N Le te. NS OO IA 3 SM Yq eo ome US OR EEE ba id LE vahy E > 2 AEN Ad e pur ce “tetas 
' Cr OUR a te wen ys Wish tra O À CRUE AL RER A (RPS BOE Teg e An Ho 5 Bes Seen tetas eri or de P rarae de Le rs a Rages 1 KA dl 
0 E a anid h E dl H Hg y GE es re Ag Rea tee ti > A os ST ro Re A. RCA LOTO oh) LL CETTE ESS DL SOLE dec ireland rodear 
G y AS da,» EN cd A A rm AT METIER TOUS Rte ret yes lod aa péri tepehahe one tn a a op Le don in 
oat > : ds nn ne ds « CL RE ES ER A rép pepe nes ds RL RE pa q a pe 
E A . A mad E TETERA eet tion pd rd eek e PES pace Ar A ata do a oe Soot = 
p : A a CPR ECC IP ETES ES K sa AR dd bee A À CHU ay Li DELLE noy e daft, me ta APE, es 2. nas LEE Aa o, 
o ' 0 Pra G N 4 LE A * «AN A A E ES A IS E] E CA a epee gered O Rott pd oe ee ee agra ea Er er En TONE A 
A 5 CPE d O ROS 7 as har» ri to e SEA LCL CORTE Sra Malas ar ma me? Ea E ee ne 
A . s , 1. A O CEN CORRE AE” one LES) HE LOL EEE" CEE CEST EN rép dd pe 2 Le pas J BRP EEE": Eee ROLE ra ied dro 
G y a Ed is i a H 4 i pl ie + D g 4 aS a ke D RAR ECS a AA A pero AY ere PRA Baria EELEE ee O LL RE EDEN E eee Spee de AS eq pere 
7 OT T i ear E “ss . e e. CET ‘#9 ee Sa ee APA Rs CR CL CCE Caro EN OS DU CL OR CPE PTS Pin. LD TERRES A APE et Et eos LR AR TROP AE nes 
a A a 0 4 Det. li DEP" TR AT epee ar) ELLE EE LR TP ver Bi Rs mh aie a he 
pal RE PAL ; : rd IO A A ON E “rey peg a pb Lees rs Bh ee eee Hr o PRA gli Cp it aa e E 
sus A AA puts O AA ES A O A A EA Y) pe bod Lela aay A ER E es pi ne le bd ree es e AO Ma mella Le ee ras 
A A i IS Y y do ra ‘ an fa ue gt ale 1 “eee A NT NS G] de tea a e aipe Ve > y E tt : E 
IS | A se AC Ha sn ` Ci ES r E E Se vus ie LA SO (moss tas sith pu proa GA o» Ph A Met yl AO uN rv Real SAR $ RE NN Der Part Sept aio fi q Ara AL an 
S o LE I. A A IS A ry AAA Fy + is % . CRETE Sjn A o A g é “ae g E ph 
E Nao ase Samad Pook REL DOUALA NA TUL LNG ABS SA E RON We COMER E E ANRE TE oe 
A CR NE G . Jen g . ` FA o g ' A AT OR ry ` > + | ny bs he ss oY EL D As s 4 Fais na hd ES AR RES SA a aan te CORRE DE en Re De CEE n e ha ane a py Corse 
| Š $ ee she Si ~ E Le i j à CL a zi i a N a OA EC + Sri tre SO RER 3 CPi edhe gr ya S yes pr Bed ÉLUS pe A eS a TAO pas w lrei 
E 0 CDR LOS ri . o UE sé ss ss» 4 E + . ETE y MA Les Ad QE LH NET cmd h ca Rese reps x et eek: RS SS o ee Ae 
"E ea E E A A RS NI gen a RTE ASC sry, CaN er ar 
K ` d "i L CENT A E a, LAS ee - a a 
do i aes i A E à. AS A CCA dl A 3 Y: : RUES MT tan 0 os ata LPS a O ds EL FENETRE a rer Mn re 
` owt 4 y dde TALA E A Ah a 1,0. PA ra [2 Pr ee e mami. se Muda Ap id mane e ho do A e Eu DR nm. CL PRES 
. Su LC AE LE A O CE Ars CT RARE D ws es ge EEEE N EAE PRE ATP En na q rs Né DE Tr o pars te 
z = É ra ' em]. . az y D 1] a. . {à 
n EM J t ts fat op CEE KE 434 ian) A Eu Re nl A Ate Stl ad Dy OS A A GEN ao UL PSE ias mia ARLES ha LE PC eee ce? ree Sree aps da 
e T E A A UE A EE RLA E EN E RAA AAE A PT ra CE ree pe iy eae ST rere ork: 
A G a QUE TL are e. Eo AL ass A TER EMTEC RSS SANT a Le COR TELE ER RE her A LP E UA o See fe 
n ó y ecm wri 0 eek i A CRT O RTS T KT vx IST Toti od ATEN QT ch att bo Gc raté Na Re VC Pr er rere Soke AE 
E CCE Meet A E A a OO MO AOS PAL “nie ae eri ee Re HE pri ye pp # à. PE ns Fer ot RTL à = 
A 4 E eg ." o L tag as CNT A OE RN ES Ia « A su wa alts Are e po ra SAT PRESSE A tp ake Re PRE à ete TT A rra pres) Ie 
G 0 e h í y e CE PEO fi sn. A E, ayi es Ei e Rs x pl i D] FRET, APO fai n° LÉ reine tien. LIT th CRE cp hat LES LA = =. 
ET A ae a # t koa rS ` A IS G - vi por sy aia Hi Li ADR ARS na ” ; b He Hick 
1” y G > d I CCE ? O E ea yinlurná midi > Un e bo Pea 
A A Chip hak] A 0 A n CCE t FLE TA 5 1 CH > ES del Sie hae jap à if WA AA ten dre LES Lor ei | PA 
i rat ies an : e RC | X EE cs LS s x - K Re YE oe 6 ae a Oe do ' pd Me to, aed a t i Re E DEIA E sas mar ni e Sos cs Sh 
E G da q LE y Sil © +» a t on g OL = + teats CRE ET EEE CIS aqi. 2am silly. éd mu. yn AO 3, e LES EN CRC EEE ie ra Aal 
y . . CA 5 y y RAR TE LTÉE AE AE LESC LT ITR) SA Des ORNE DLL TE Lb ART er e di EE 
E ae. E . + : py A OS TO Tr HU CORR NA Re A PP pr el ARTE DATES A Ca: id a PTT ORNE Te or fees TE 
E 7 p . O . . «e y CR LE CE ry 7 A y à “J 
EL G units ee E AN A Re a AA Aca eRe wh Saba TREY LS CES E A ra Er dd E 
. . o a A . . g ’ . SS TL * Th > y ES q oars re CLP, TE a £ 
a f J o y a Ne eee mr E È e E ad e EF i3 H Era fy FR ue CLR Pe ee Men RIT y Dc de is SNS Lima! Ge bert a A ER sax reg 
O 1] g 7 a ES TT LEA | we Ne CRC ws ad i. S CLR bo A À x Fi ch f 5 h per 
' h ar È 7 A y Pa A soo vi een af e a? re A A aS ry A ONE tet fa DETTES en DATE TREA de AH “3 SAT Ae OA LE E Sie PART nr ‘= 
A s 0 E Pa (E E E) oo . 5 TE A A pw ow, Ae iat) bd Weer oss AS JO e ES ds de O T DE . A LR wiy o FA pe > a : + CARPE L 7 
E QE l d LR a a ache CC LION, PAC CeCe RCE DR RTE TS stef CCR AAA LAN CR NE SH AT y de SA =) o n a 
4 ' 1 g A st * ré t A , 7 , LE E r H 0 Ps Ta E NS pe € 
F a P T E At A ae Sie RE hr CORRE ES AREAS Pie 
. E A ' g A x A a A ki s 4 ty e ET teat “' y ” ye 2 fAs ti Arda Sate AC aan 0 i ah: The e a 
rae Parl . ho» ss PC A ET ANO F PNEU CA > DE E TC A N el A Re 2 A re À zag ie) per PA eed E CAES 
MERE H y Ot: A OL AS] Ps ha oa RE ty A ARH e ES A a y patente Lae 
X oe ia y À i y 4 oder ees rea Jet Ma Re MT nue (ie CH REX da ce y the e Ne A ee PR O Tr CS Con 
: ` SL E S TA CE Mess a A IO AS ts E 7 ARE A a AA ae QC PART rit 3% Le ty ed e A 
p : ... AS E A E beset Os Sant E SA ae Fi KT CUA pioneers 
i i RS , re a ; z 4 3 Le Ton e 4 DS ne Eco E eee! La ¿y e, = DB ir es” ejes 
En Y E AS A OS EN ar A E AOL 24 9, ie EY A a Rd a ates EEE CH rer PA 
SR o - fe RO RE TN EE IR E 
Sia : A n°? A À fn 5 ce 5 T 
= E i E a E + a AOL DRAC LEA si > yi Paan anita? A SA ds O, Ter P A f, RA mr 
A 2 Bed ù RU A 4 He Va At LSS CARE EE a e ni A ee ee 
A e 3 Meare E ete NC tit PA RL A F PRO aye PAU KE «POS SA E PRE e 
A le DAC 5 ES aa y ES at RATAS E TOUTE A e 
A 7 A 5 g 5 enn. DL Pray O p CR Pe! a r < sy 
2 Si à 1 ina J nd F š di Pe igen re Pre eae OR $ net $ A F Re ue ? ES bin 
a A ae A E a ai fe et} RUE E Le U = Ce e RU MEL 
E el ba PRE i ente CR ~ SE 3 LPC 
d A a j a E E ak = LE At uin ITA TEE TE rie es ofS La A AS ae Fy 
we E 0 su OPEN CPL DE DCE e RCE RE +a : COPAINS UTP PES REINE NC LE UA ora Ms HOUR LR ete 
“Ty 0 G G 4 de Soy lata Sy VTL lil OG vER RER er HSA T ES es rie e ee 
Ji A 0 0 A | H e io zo e a? EA A TA ES en fafa RE TRE S entes! + tah PS to Y ea ES O 
Ge i A , P . LN A LY A z e e o: FANS CURE PTE CORRE CRAN PE CET r 5 Y” FE Er, = ng ie re es 
E k CE Y i Wa, Pao . A TO E NARRA D E sare af than 
vs i Coe enh” E oe 5 ds ME ds j Ae A RAS CONTE Cr sor E de 
a vue CA 8 a y y IR ES O oe j AE K EC PRET PC Er pee a ES A AA 
P A a 7 ' ' UN] a ue AA E P Ur re a Ni Ta LE y » LEE RAYON Beas = AP Se PTT Dr PT ET 
iig 7 g E Wie! ie . Leen | ’ ar LL NE ` POS io po ES ur? SALE ett | site AS A 
A CE E ANA Po PY COURANTE ON Mes RO du t i k a $ O 
g 0 13 , ES 4 EERS P CH LT 2.2 A é F 
. . . A i = °F n . A A y n 1" a » pe Shas Py P mi EF-S ER s P 
Û j 1 be i 3 ee a ds AT Pr CARE 
A A A Ae che à i s a f es O A e O mt g g, is LPS TR È er 
. y is p D a Pan ae 
5 5 g A Pry reat) E 3! CRC Seale cry ee ey sae WE PRE 
a n LU 5 .. ¿. 3 ef, A J l 3 U eiT) T y y LEEN AAA 
= A | TE Ar g * CT £ ` ir CERN, nd" ERA A CES % 
$ A te A oP E ee uae Pee OCR CEE ET re ty IA Ad: 
A D A Sr a H 5 
o g AP k dr LN i a Al y 59 .* LE RONA Ne J 1,3 O CM a rau! * t, tt s tur Ln ss 
inet A ... 7 a CRC i E ate ca. AL ey i der as a PRINCE à 4 es A hla Eh 
po oc ee HE A Ra td a3 A A ni De ES CORTE UE A PO q a Pa 
A ei . A eee H 0 y E E g Pek E L U, { As EN ME € = ow ey Ob Yee te Poe OS Scie 
nepeg: g P à CE CN] FE he SP ti o ee ac “AE ai ón š N 4 a Steaua ld... AOL AS, as] A < 
G é ) CE Aoa ya et je e “rt i s RC sae ea 4 TRD SA FERET Kae Lu se “A 4 a É y: 
RS ig A Mehr M 5 9 A LS ER A 7 AAN PE 
y y i O Y tos Set ES D ate Ne EL HA PR A Bao Ay Me Mee Gre ty ge ake LP re ie A E Lave 
A A y CE A te ng ' E CR Lee A NS y CP ET] OR CAE Hs PAL MERE dr PEE ee dd cp alg 
' à 4 0 o E O ‘ E d 6 CR ETES pr P FD VON I NE TE eu g ne rag LE oe ee 
y g A AO J bai hal Tae} COR CORRE DC ROIS ECT a ets FE" 2 ved v Sos H ane PP) 
A A 5 g AO E AS LE LCR ‘9 fi F DE TOR À E RS ES Er TA ne cd pan LA Ur LE RES pr > 
e eer | ete + @ ' CON 0 ry AU 7 CE AT 7 2 .. ia or te DD LOL a 
i Bl 5 PES H Fa als RUE 0 Vy LR LE a j F. 4 by RARER TE dus LA À ops are ES à CA aa LL pre E 
i i o o x 3 ye Hai 2 i $ EE AS E Ey cae Eu er E te aks y AUS VE CE A nine , 
E ig CE DURS CR TE a DE EC DT 0 A E f Pe sF By ae Sar VE T HAR ah AT po Po RL AE ropes 
D ed ' y A NC GA A OR edo 7 ar E RE PEA AS $ AA co PEPITO DAS 
A as" À r 
i A IL ROC y A Aas ur QU HE 0.) carne ty re tor QUO as AL RTL ST CS PE EE RO PT ERA 
i n CAE RCE ‘y e F .« re y E E 1. Han À UN 5 A 2.4 arte LP perds TA ROE O TT ALORS pops 
à G n i se re cia ere i pli! a. i D SC de PT ART Cr) fs REE At ele RE ae a a PET DT A Ary y 
$ n LR RO CN pore? r, LP AT ANRT NICE A DRE ae LE ph SUR COS CET CORP EE EN EPS re beara 
i E LU Pr 4 hos > nr Fete! OS US be ue a LES 
a n ATA | sy a : y la P Qu FA y E = SO rt Py CPR PE CET Bene LAE EO ees rae 
ag , IA pi A 4 r REE E T MARS T Res noe RES PA ART GES Poll teed 
On cane uy CREME Jé n Ex A IA EE PL TI Ts sf EEE ST es 
E T) te # Zn CAES ma PE pe sf panas Rte RARE RUES 
IED ae E = g HERTZ mb LA WE e fe Pe 
4 fr LE: SE RS oad) 
CPL CLA De Tate PRE a ES ELA As PPT MR RES fl dE pjs 
ciara je DRAC A ONE Rete rae yi RO he CRT ETA PP ET AP er we 
AR O A A : CRIE ST va A eb. PA e A ER A dos 
O IRC BO a eo MC ARR NT PS | PNR UN NT PER RE ER IE e 
K aks nl M aly? wap ta FLE Fra ns AA CON ss Pe Dn de Ki nn Beto mt ds PE A Gers LA A | POUPEE PE PTT 
i e ae e a Fe z nee i” CEE te ee ET k f PJ nes Ada RARE LIEN RE LEUR 
COR AC | CRUE pts di. ARA TRS 
6 z. AAA | 


s 
DR LE DA AUX 
CLAIRES 





4 R m ng PES CNE SL 
(YA or nee LEA rr Aires c= ATA 











GO fps gia Bei 
y 49 A Wor | 40% dr CARD LISE À y Fe PAT Var PRET 4 K 
ese tate CE HS o r} PEA LUE CU P EE Cr mie EI: A Pir eters pple die rn Cg ag 
ee ECS | Ait D Y AA A 1. OA OO IRSA UE ES thd HS her ses f 
OT E EEE e / ne D LIA AIDA CA rair Pe ase Ao TS aig: NE LUN AS ET HE pees LE 
y Ve. CE) ne, apa E O P R Of Ea RAS PEER ES oe 
CEA +. CR A TE ERT E EROP ee SCR ee ee pe HAVE COCA s 
a AS E OS E 


E BS FE OIA 8 





br gr i ed Be ss Raden at toe) Tes 
FRAPPE EUX CCE E OT AUTRE y Ae 


LA IEA A 
$. Pide AN 















on 
SR AS y TIERNA RS PORTÉE Ve 
IA U CREER de A E ne fi rit Pg oe) il hina BLO at le LAA AO 
O Aé ES is Cr E E oe EE EAP LR Sorts uT ate RC ERA te RTE CL ao PPA ere Pree 
cn E E oP NV UA Det itn AT hat ee it ES vee x PTE PRET EEE ’ 
Py #, a AER IRE A AA ‘1: Ve ene We ree? CPE EL CET EN pai As P 
F s BR ‘ rary Mn AN AS LE NT TER à 7 RER D ER AS A di ia a aa h 
D MORTENS ROLE DER es ES SA DSL e IA fa e AO PRA a 
Z re RAD NN f "Vat the Re DE CLP DD se LON GS api ee fA 
A EF COTE Y A A RA EI en 
e o Li 
dto o. 5 





CRUE AL 
oie eu € RUE] RO 










RAA GE 0 
Let AS "Ad Air Ps pr 
g Y $ p ARE ak ee) A TT E ER N rt 
4Y A E A a T 


LI 
a yo PEEN) Sa HAE EP EE ZT Mer 
TAO 14e. Dan Pie he Lge EA, RANCE 










































































































LUE AAN A 
| BLOT LN De ere tah aan ARR EAP NS AE ad RAR 
4 y o M OL A e PER TIRA A PEN $ y COE T p CASET aw ey E ae y 
0 AE ne PL Er CE te EU AT ta +, TA Nr PAID dee KE RSS REY ET SE 
ro À Ft 7 gi A AE w ee fae Fe fi KT CU ae ye 0 E JS ASC ML CRIS here ne othe tato 
a . y oars ev v Fe À wae 0 f pay 
ies 07 k w PI TN EN y DD A pe N AS O ie LS Cat PE arr AE À RES 
A MVE AE CUT UP Pe Par RRE e 3 ae y = i jee PEA a AAA 
E EE OS E RI ODO CA pe CAE A me pein P 
NN LES ET SA ER OT PE Ot LENA Ee 
i z 3 s A 4 4 G x L i k e =! #7 a 
NX cda! CRU Fo CEA PO iy ye AE ne CORSA CE OCT RTS > y PE ns Ae RE LE PASE ia A | EE és PAS 
aa a RO AS E ES ELN ULSA AE AR ee te ae 
PT LPS RAA : UC PAT E CRE SUR DO LE CE PRE PS PTS Faty O EEAS TARE Pi EP A T 
e 0 IES E TCAE PA £ oe bì ET La pee PE vtr wy? COLA TE E OT . be iE ara wh # A SAA He. en = feta 
, A Pid JA CA D Fe LAS ne CRT: . O AR SY AI IDA we % Peps thet Ga PU fs AS Eh 
a oF Ri R E E E r ASE Ae w ARE EST RADIO EN Y E E Ae PONS ake Oe sn LIT pips ral 
ALT EU E UE Er RE EAS NAAA eee MAS A he ELU ver a 
EEEE, k KA NEUTRE CNE x DUR CET PET k t a pe RMS RARES en PE Re re T 
sie PRES RC ON" q PAP PE e. ALLEZ OT TEE i 
IN A EI TO O rs LAN) PS 4.4 Py CA LRU 
ta ba Demme E T P P , r A = Arie Cp Nr a A eS EA ia > 
He A Ae et dd NRA ES - ge eo MCE sly prado A E AS CARNET CADA 
MANE O 1 nose bs OA AS TE E ie be we TS 3 DER bad ay Le y e RE TE AAA od C2 
CPC Hi G f A i E A ng AR hed a: 1 VIE E TT LARA ee 
O F à AA p S Dos te re athe ets QE PEM Tee Ce Pon aa ie ‘was Peek tre er rae Pad tah Seedy bd 
i ard LE j 7 f rT a es See ee ee RT eo See A RES ASA ALA Ce Fs e al ere gle, 
ap rr’ > rid r Py 4, N Pi oom E eS OA FA EA A Pid E ie pl Aa ee Srey fe A as CE CAR 
, K TL EN AA Le x Cf EPA Fi + 5 aa 
0% ae y EU A LO E E Ra. 19 AT REA A fa AN IRAN ASEOS A EE CLP 
PUR RCE DORE CE RU MITA. COCO ES pu TT CORRE ATEN UE DE ERA A AAA ay” e ds ie a 
AA A AAA DAT EAS E UE ge Seo ON DE PR AS ey es es a y TS LAA 
' 1) . n ry ry Es TE PT BN es 4 4 7 
vi 7 ts ye EC EU as PN i. Ey > ur ik x A nt ERREA LA AW cy re RAA y 7 +, ARA ACER ee dE Foal oh N. 
J CR | EET) 3 E 3 
SARA ETC SORTS En LR E TED SPA aS AA PNR ESSE CR ee 
a A ie 0 ce Ae 0 ba LAO PARAS ni É a Ce a AA ey? r NES EE lan E Y NM AO A: PH 
TA A OL TA REC à PA ASS TR ARO AAN ALO EPA A a a 
EE FEE a y A AI SE E E E RER AER AAA A RAPT SS ve f 
of id COLLE RUE OS er Te oe A TE W DS i, QUES EER RETE ES Al er A ES 
Nr La CC ERA MARAEA 0 TARER , i: ARR KOIR Hee os AECA ST E 
F Ë ES ae de Ca F ra OO On per A f be 
es HS e ne, ; ced 3 HE + rp ESAS PE q AR (y DOV i ego ANSE EST CLR EAU se + s kel, Lig Sb 
LE Ars HE PEEL Me EOL E CIA Pa Ren tee E ET be RAILS ee 4 pe dit 
AS e SO rs do Hepat wh ade Tw, GES AE LR A LE RENE ES En Li E 
a aL a Mare) oT lek r O E II A ANA EA TNA IA M 
fore PE LU PE a es if AS Fr CCOO ce he Ls ¿pr À 0 ea pl haba o EPN 
y 4. CA ER CR PE B P LA 2 
hae a coe E E COUNTER LR ea A Ÿ. OS ets rAd Ed A AN 
1 ent Lo vs e wns a EAU APCE, LR yt AO pA q LA SAR q o A EA AAA hth 
eT eg ered E E ll a y e do nds 
q + PNA ate Per CRUE REX ECO A A A ty mE AS AO h DA SR 
iy Ade EN PURES A E SCORE PTS ASAS BER AAA NA 
A 0 be a = A ACA ' O it ? jasi 5 AS A E COS NN NT CS ARR a 
A { ; Dae g A AR y G a AU Tor CC fps A ne DE er AE: PE E QU PACE vé De E IAS 
Ea 0 5 OA aes | . g AN ” ys TT EE Y x VI 
4 + ’ RE rhs acer PR ro s rh KAN AX LE Re IAN en ernie SS ee) eo DON AA 
ma 0 A 3" st A bos PECE; CL SR E 2 = Wages aA e UE A e PA id parado “Re 
A EF) A E PET bts, CCS Les RL NE ON AE a art ane p x A E Ded 
à G A A E i ER rag ae Ler a a w a > A O ñ ee re EG m A o E A RATO, 
‘ 8 g ray | , A] P) E . ch 2 + 
g 0 1 n 12 1 E A HT, Le A Tu, 
LI ` 




















NPS-OC-91-00 


NAVAL POSTGRADUATE SCHOOL 
Monterey , California 








THESIS 


EVALUATION OF ABSOLUTE POSITIONING US- 
ING THE 
DEFENSE MAPPING AGENCY’S GASP PROGRAM.: 
by 
Dennis Bredthauer 


September 199] 


Thesis Advisor J.R. Clynch 





Approved for public release; distribution is unlimited. 


Prepared for: 
Naval Postgraduate School 
Monterey, California 93943-5000 


1257807 





NAVAL POSTGRADUATE SCHOOL 
Monterey, California 


Rear Admiral R. W. West, Jr. H. Shull 
Superintendent Provost 


This thesis is prepared in conjunction with research sponsored and funded by the Naval Postgraduate 
School (NPS). 


Reproduction of all or part of this report is authorized. 





Jnclassified 


ecurity classification of this page 


REPORT DOCUMENTATION PAGE 


la Report Security Classification Unclassified 1b Restrictive Markings 


2a Security Classification Authority 3 Distribution;Availability of Report 
2b Declassification; Downgrading Schedule Approved for public release; distribution 1S unlimited. 


4 Performing Organization Report Number(s) 5 Monitoring Organization Report Number(s) 


6a Name of Performing Organization 6b Office Symbol 7a Name of Monitoring Organization 


Naval Postgraduate School (if applicable) 35 Naval Postgraduate School 
Monterey, CA 93943-5000 Monterey, CA 93943-5000 
(if applicable) OMEN 


Sansa nd ZIP code) 
11 Title (include security classification) EVALUATION OF ABSOLUTE POSITIONING USING THE DEFENSE MAPPING 
AGENCY’S GASP PROGRAM 


12 Personal Author(s) Dennis Bredthauer 


13a Type of Report 13b Time Covered 14 Date of Report (year, month, day) 15 Page Count 
Master’s Thesis From To September 1991 161 


16 Supplementary Notation The views expressed in this thesis are those of the author and do not reflect the official policy or po- 
sition of the Department of Defense or the U.S. Government. 






17 Cosati Codes 18 Subject Terms (continue on reverse lf necessary and Identify by block number) 
Field Subgroup | GPS, GASP, Absolute Positioning 
eae LLC 





aoe | OO 
19 Abstract (continue on reverse if necessary and identify by block number) 

The Geodetic Absolute Sequential Positioning (GASP) program, as utilized by the Defense Mapping Agency (DMA), 

processes static GPS measurements collected with the TI 4100 GPS receiver to estimate geodetic point (absolute) positions. 
In this thesis, the GASP program is modified to accept data from different receiver types, the estimated point positions are 
compared to positions produced by the Transit Doppler positioning system, the between-receiver estimates are compared, the 
erence between estimates using the broadcast and the precise are examined, and the effects of Selective Availability as- 
sessed. 
During the Monterey Bay Precision Positioning Experiment (MBPPE), conducted in the Winter of 1990-91, a large set 
of static GPS positioning data was collected with four types of GPS receiver; the TI 4100, the Trimble 4000ST, the Ashtech 
LD XII, and the Magnavox MX4200. Additional static GPS measurements were obtained with the TI 4100 receiver at a 
‘eference site established to support the expenment. A third data set was collected after activation of Selective Availability. 
Measurements collected with the TI 4100, Tnmble, and Ashtech receivers were subsequently processed with GASP using 
broadcast and precise ephemerides to produce point position estimates. In order for GASP to accept the data from the 
Ashtech and Trimble receivers, the program had to be modified. 

The positioning results obtained are analyzed for accuracy and precision. The accuracy of the GASP GPS estimates is 
letermined by comparison to independent estimates obtained by the Transit Doppler positioning system. Precision or re- 
yeatability (1.e., consistency of the estimated positions) is also examined. 

Analysis of the accuracy and repeatability reveals little difference between the positions computed for the three receivers 
ising the precise ephemeris and that all three provide good agreement to the Transit Doppler positions. All three receivers 
re capable of providing geodetic-quality point positions. It is also clearly demonstrated that the precise ephemeris does 
| roduce a more accurate, higher precision solution than the broadcast ephemeris. The activation of Selective Availability has 
ubstantially degraded position solutions available from the broadcast ephemeris. 


0 Distribution/A vailability of Abstract 21 Abstract Security Classification 

À unclassified/unlimited O same as report O DTIC users Unclassified 

3 Name of Responsible Individual 22b Telephone (include Area code) 22c Office Symbol 
,R. Clynch (408) 646-3268 OC/Cl 


1D FORM 1473,84 MAR 83 APR edition may be used until exhausted security classification of this page 
| All other editions are obsolete 









Unclassified 


p= © 


Approved for public release; distribution is unlimited. 


EVALUATION OF ABSOLUTE POSITIONING USING THE 
DEFENSE MAPPING AGENCY’S GASP PROGRAM 


by 
Dennis Bredthauer 
Civilian, Defense Mapping Agency 


B.S., Worcester State College, 1982 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN HYDROGRAPHIC SCIENCE 
from the 


NAVAL POSTGRADUATE SCHOOL 
September 1991 


A 


ABSTRACT 


The Geodetic Absolute Sequential Positioning (GASP) program, as utilized by the 
Defense Mapping Agency (DMA), processes static GPS measurements collected with 
the TI 4100 GPS receiver to estimate geodetic point (absolute) positions. In this thesis, 
the GASP program is modified to accept data from different receiver types, the estimated 
point positions are compared to positions produced by the Transit Doppler positioning 
system, the between-receiver estimates are compared, the difference between estimates 
using the broadcast and the precise are examined, and the effects of Selective Avail- 
ability assessed. 

During the Monterey Bay Precision Positioning Experiment (MBPPE), conducted 
in the Winter of 1990-91, a large set of static GPS positioning data was collected with 
four types of GPS receiver; the TI 4100, the Trimble 4000ST, the Ashtech LD XII, and 
the Magnavox MX4200. Additional static GPS measurements were obtained with the 
TI 4100 receiver at a reference site established to support the experiment. A third data 
set was collected after activation of Selective Availability. Measurements collected with 
the TI 4100, Trimble, and Ashtech receivers were subsequently processed with GASP 
using broadcast and precise ephemerides to produce point position estimates. In order 
for GASP to accept the data from the Ashtech and Trimble receivers, the program had 
to be modified. 

The positioning results obtained are analyzed for accuracy and precision. The ac- 
curacy of the GASP GPS estimates is determined by comparison to independent esti- 
mates obtained by the Transit Doppler positioning system. Precision or repeatability 
(i.e., consistency of the estimated positions) is also examined. 

Analysis of the accuracy and repeatability reveals little difference between the posi- 
tions computed for the three receivers using the precise ephemeris and that all three 
provide good agreement to the Transit Doppler positions. All three receivers are capable 
of providing geodetic-quality point positions. It is also clearly demonstrated that the 
precise ephemeris does produce a more accurate, higher precision solution than the 
broadcast ephemeris. The activation of Selective Availability has substantially degraded 
position solutions available from the broadcast ephemeris. 


THESIS DISCLAIMER 


The reader is cautioned that computer programs developed in this research may not 
have been exercised for all cases of interest. While every effort has been made, within 
the time available, to ensure that the programs are free of computational and logic er- 
rors, they cannot be considered validated. Any application of these programs without 
additional verification is at the risk of the user. 
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I. INTRODUCTION 


A. THESIS DESCRIPTION AND OBJECTIVES 
1. Monterey Bay Precision Positioning Experiment Overview 

The Monterey Bay Precision Positioning Experiment (MBPPE) was designed to 
assess the performance of some commercially available Global Positioning System 
(GPS) receivers and the processing software developed to support the receivers. In 
particular, the static and dynamic positioning solutions of the Trimble 4O00ST, Ashtech 
LD XII, Magnavox MX4200, and Texas Instruments TI 4100 receivers were examined. 
All but the MX4200 are geodetic quality receivers. The Trimble, Ashtech and 
Magnavox models exhibit some of the latest developments currently available in GPS 
receiver technology. The TI 4100 is an older model but is still widely used throughout 
the geodetic community. 

The overall objective of the experiment was to acquire a large data set of GPS 
measurements in both static and kinematic modes with the different receivers, to process 
the data by utilizing processing techniques and software relevant to a given application, 
to interpret the results from a position accuracy and error analysis perspective, and to 
evalute the receivers and processing techniques based on these results. 

2. Thesis Overview 

The purpose of this investigation is to determine whether meter level point po- 
sition accuracy is attainable with the TI 4100, Trimble, and Ashtech receivers. The 
Magnavox receiver is not considered in this study because it is a single frequency re- 
ceiver, whereas the other receivers operate on dual frequencies. This is a serious limita- 
tion in determination of point positions since the dual frequency correction cannot be 
applied to correct for ionospheric refraction. 

This study is concerned with one aspect of geodetic positioning, the absolute 
determination of a point position in some commonly used reference frame. There are 
two general categories of positioning, static and kinematic. As may be inferred, static 
positioning involves determining the positions of stationary objects and kinematic posi- 
tioning the positions of moving objects. Static positioning may be divided into relative 
and absolute positioning subclasses. Relative positioning is described as the estimation 
of the vector (baseline) connecting a known station to an unknown station. Absolute 


or point positioning may be described as the estimation of the vector connecting the 


origin of a global reference frame with an unknown station [Ref. I: p. 1]. It then follows 
that relative kinematic positioning may be described as the estimation of a vector con- 
necting a known static station to an unknown moving object. Station positions deter- 
mined by absolute positioning methods may be used as mapping control points, datum 
transformation ties or in other mapping related applications. Once the absolute position 
of a station is determined it may be subsequently used as a known station in a relative 
positioning survey to establish the positions of additional stations. 

The tool used to produce the absolute position solutions was the Geodetic Ab- 
solute Sequential Positioning (GASP) System of programs developed by the Defense 
Mapping Agency (DMA) [Ref. 1: p. 1, Ref. 2: p. 487]. GASP estimates the absolute 
position of a point in the geocentric, earth-fixed, cartesian, Conventional Terrestrial 
(CT) coordinate system in the World Geodetic System 1984 (WGS 84) reference frame 
(See Figure 1). One axis of the CT system passes through the intersection of the 
Greenwich meridian and the equatorial plane. The third axis passes through the Con- 
ventional International Origin (CIO) which is the average position of the earth's rota- 
tional pole for the years 1900 to 1905. The second axis is orthogonal to the first and 
third axes in a right-hand sense. [Ref. 3: Sec.5.6] 
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THE CONVENTIONAL TERRESTRIAL COORDINATE SYSTEM 


Figure 1. The Conventional Terrestrial Coordinate System 


DMA has been processing GPS data collected with the TI 4100 with GASP for 
the last few years and it has demonstrated that geodetic quality point positions are rou- 
tinely achievable. This means that a position solution with a standard deviation on each 
component of less than one meter can be expected, usually with data collected over a 
span of about four hours at a data collection rate of 30 seconds. Because accuracy re- 
sults using the TI 4100 data with GASP are well documented, a benchmark had already 
been established that would provide a reference for the results we would obtain. [Ref. 
L pS] 

a. Thesis Objectives 

Before this project, GASP had the ability to process data collected with the 
TI 4100 receiver, using the Floating point, Integer, Character, ASCII (or FICA) format 
as the input format. One of the goals of this thesis was to modify GASP such that the 
programs would accept the measurements of any GPS receiver type. The Receiver 
INdependent EXchange (RINEX) format was chosen as a second data format through 
which receiver measurements could be entered into the program [Ref. 4]. The RINEX 
format was selected because it is one of the more widely recognized and accepted of the 
exchange formats currently in use and because both the Trimble and Ashtech data 
processing packages have programs that will convert their data to the RINEX format. 
It was also desired that the GASP programs retain the capability to process data in the 
FICA format for which GASP was originally written. 

In addition to modifying the GASP program to accept the data from the 
Trimble and Ashtech receivers other objectives of this thesis are: 


e To compare GPS position solutions produced for the TI 4100, Trimble 4000ST, 
and Ashtech LD XII GPS receivers to the solutions obtained from an independent 
method.l 


e To compare the receiver position solutions for the broadcast ephemeris versus the 
precise ephemeris 


e To examine repeatability (i.e., consistency) of results 


e To examine the effects of Selective Availability on the position solution when using 
the broadcast ephemeris 


1 The Transit Doppler positioning system was the independent method used in this compar- 
ison and the positions derived from this method are taken as the “true” positions. 





B. BACKGROUND 

1. The Global Positioning System (GPS) - Fundamental Concepts 

There are a number of satellite based positioning systems being used to establish 
the position of an observer on or near the surface of the earth. The Global Positioning 
System (GPS) is one such system [Ref. 3: Sec. 3.0]. The NAVSTAR Global Positioning 
System is a passive navigation and satellite positioning system operated by the Depart- 
ment of Defense (DoD). It was developed to afford the user instantaneous three- 
dimensional position information anywhere in the world. GPS exploits simultaneously 
received radio frequency signals to determine range measurements between satellites and 
earth based receivers. These measurements along with a knowledge of the satellite posi- 
tions can be used to solve for the receiver coordinates. 

Through recent advances and upgrades in hardware, software and the develop- 
ment of new techniques in data processing, the systems applications have expanded 
dramatically. Particularly in the areas of high precision surveying and crustal deforma- 
tion studies, GPS has rapidly become competitive with other positioning systems and 
techniques [Ref. 5]. 

The strength of the GPS system in geodetic work lies in its relative affordability, 
portabilty, ease of operation, and high accuracy in comparison to other high precision 
positioning systems [Ref. 6]. For example, the Transit Doppler satellite based position- 
ing system requires two to four days of data collection to produce a geodetic quality 
point position. GPS, on the other hand, is able to produce a geodetic quality point po- 
sition in only four hours of data collection [Ref. 1: p. 5]. From a mapping and charting 
perspective, the fact that it provides position determinations in a unified coordinate 
system may be its most important feature. 

a. GPS System Components 
GPS consists of three primary segments: the satellites, the ground control, 
and the users. 

(1) The Satellite Segment. The GPS satellite constellation, as of mid 
1991, consists of 16 satellites deployed in high earth orbit (about 20000 km altitude) 
configured such that a minimum of four satellites are visible to the user at a given time 
with from five to seven satellites typically available. The satellites are arranged in six 
orbital planes and have an orbital period of about 12 hours. It is proposed that 24 sat- 
ellites eventually be deployed to ensure that a minimum of six satellites will be visible to 


the observer. 


Each satellite transmits two radio frequency signals in the L-band, 
one at a frequency of 1575 MHz (termed L1) the other at 1227 MHz (L2). The LI fre- 
quency is modulated by the Coarse/Acquisition (C/A) code, the Precise (P) code and the 
Navigation Message. The L2 frequency is modulated only by the P-code and the Navi- 
gation Message. 

The C/A and P-codes are binary, pseudo-random noise (PRN) codes. 
Pseudo-random noise codes resemble true random noise except that random noise car- 
ries no information whereas PRN codes are generated to be predictable and carry in- 
formation. In this case, the information is used to determine the time of signal transit 
from satellite to receiver. Identical codes are generated by both the satellites and the 
receiver and cross-correlations (determination of the scalar product of the code sequence 
with a time delayed copy of itself) of the incoming satellite codes with the receiver gen- 
erated replicas are performed to determine the time reading of the transmitter clock. 

The C/A and P-code modulations shift the phase of the carrier L- 
band frequencies by 180 degrees. They are essentially a sequence of positive and negative 
ones superimposed onto the carriers at frequencies of 1.023 MHz and 10.23 MHz re- 
spectively. If the code value is minus one the carrier phase is shifted, if plus one there 
is no effect on the carrier signal (See Figure 2). 

The C/A-code sequence repeats every millisecond yielding a C/A code 
cycle of 300 kilometers. There are 2'° C/A code chips in one code cycle so that the length 
of one code chip corresponds to about 300 meters. Since it is possible for receivers to 
measure fractions of chips, the C/A code may be used as a medium accuracy navigation 
signal. C/A codes are exclusive to a particular satellite making it possible to distinguish 
between signals received simultaneously from the satellites. 

The P-code sequence repeats every 267 days and is subdivided into 
38 seven day segments. Each satellite is assigned a one week segment of the code. Thus 
all satellites can transmit on the same frequency and still be distinguished from each 
other. This weekly subdivision creates an identification system based on the PRN seg- 
ment assigned to a particular satellite. If, for example, a satellite is assigned the seventh 
weekly segment of the code sequence, it is identified as PRN 7. Codes are initialized 
once per week at Saturday midnight. Since the length of a P-code chip corresponds to 
about 30 meters, the P-code supplies a more precise measurement than the C/A-code. 

Critical satellite hardware components are the on-board atomic 
oscillators (highly stable and precise cesium and rubidium clocks) which control the 
generation of the carrier frequencies and code modulations. The signals are coherently 


Code Modulation 





Figure 2. Code Modulation 


generated from the same fundamental frequency, designated f} = 10.23 MHz. The two 
carrier frequencies are multiples of the fundamental with f,, = 154 f and /,, = 120 f, The 
C/A-code modulation has a chip rate one-tenth the fundamental frequency (f}/ 10) and 
the P-code a chip rate equal to the fundamental frequency (f). 

The Navigation Message is a low data rate message (broadcast at 50 
bits per second) that contains information on satellite health, satellite clocks and the 
Broadcast Ephemerides. The contents of the navigation message are updated hourly. 

Before the receiver position can be determined, the positions of the 
orbiting satellites must be known. The broadcast ephemerides are a set of predicted or 
extrapolated orbital parameters which define the satellites position with time and are 


used principally in real-time positioning. Satellite tracking information, received by five 


ground tracking stations (identified in the next section), is used in a least squares ad- 
justment to estimate the Satellite positions. The positions of the ground tracking stations 
are held fixed in the estimation process. The broadcast ephemeris information is valid 
only over a specified time interval (about six hours), not over the entire orbit. The main 
reason the navigation message 1s updated hourly is so the broadcast ephemeris remains 
current. 

The Keplerian orbit description contained in the broadcast ephemeris 
is the means by which the conventional terrestrial coordinates of the satellites are com- 
puted. Six Keplerian elements and a reference time are needed to completely describe the 
Satellites orbit in the CT coordinate system. Five of these elements define the satellite 
orbit and the other describes the position of the satellite in the orbit as a function of 
time. The other broadcast ephemeris parameters transmitted in the Navigation Message 
describe the deviations or perturbations of the satellite motion from the smooth ellipse 
defined by the six Keplerian elements refered to above. For the complete list of the 
broadcast parameters and the computations required to convert the Keplerian orbital 
parameters to coordinates in the conventional terrestrial system, see Appendix A. 

If more accurately determined satellite orbits are required than are 
provided in the broadcast ephemeris, a post-computed or precise ephemeris may be ob- 
tained. The Defense Mapping Agency and the U.S. National Geodetic Service are the 
agencies responsible for generating and distributing this information to the user upon 
request. Because the precise ephemeris is computed after satellite observations are made 
and is an interpolation of satellite position based on a least squares adjustment using ten 
ground stations rather than five, much better accuracies for the satellite position are at- 
tainable. The precise ephemerides, as supplied by either DMA or NGS, differ from the 
broadcast ephemerides in both data presentation and in the frequency with which satel- 
lite positions are estimated. The precise ephemeris furnishes satellite positions already 
computed in the CT coordinate system and gives the satellite velocities in the directions 
of the coordinate axes. The CT satellite positions and velocities are determined every 
15 minutes in contrast to the hourly updates of the broadcast ephemerides. 

(2) The Ground Control Segment. This segment consists of five globally 
distributed monitoring stations that track the satellites and transmit tracking informa- 
tion to the master control station in Colorado Springs, Colorado. It is here that com- 
putations are performed to provide the broadcast ephemeris contained in the updated 
navigation message. The updated navigation message is then uploaded to the satellites 
for broadcast to the user. The navigation message, as previously stated, contains infor- 


mation on the health status of the satellites, and information on the satellites orbits and 


atomic clocks. The five monitoring stations are located at Diego Garcia, Ascension Is- 


land, Kwajalein, Hawai, and Colorado Springs (See Figure 3). 





Figure 3. GPS Control Stations 


(3) The User Segment. The user segment is composed of civilian and 
military users exploiting some type of GPS receiver and receiving antenna system. GPS 
receivers have the ability to receive the L1 and L2 signals from a number of satellites 
simultaneously by devoting specific channels to each satellite signal. Once the signals are 
channeled they can be individually processed. The three most commonly used types of 
GPS receiver are the code-correlating receivers, squaring channel receivers, and receivers 
that incorporate elements of both of these techniques. The code-correlating technique 


is normally employed to decode the broadcast message and also to provide pseudorange 


measurements from either the C/A or P-codes. The TI 4100 provides pseudorange 
measurements from the P-code while the Trimble and Ashtech pseudoranges are derived 
from the C/A-code. Once the carrier signal has been demodulated (i.e. the navigation 
message and PRN codes removed), the carrier frequencies may be processed to provide 
carrier phase measurements. The squaring channel technique provides carrier phase 
measurements but cannot be used to extract the broadcast message (satellite position 
and clock information must be supplied externally) nor provide pseudoranges (for details 
on the pseudorange and carrier phase measurements see the next section). 

Essentially, code-correlating receivers operate by internally generating 
a replica of the incoming code modulated satellite signal. The start of the code replica 
will be offset due to the propagation delay between satellite and receiver. The replica is 
then incrementally shifted or cross-correlated with the incoming coded signal and the 
clock controlling the generation of the replica signal is corrected to reflect the shifts. 
Once the replica is aligned to the incoming signal, it stays locked to it. At this point, the 
replica code generator clock reads identically to the satellite clock and the signal trans- 
mission time is determined by differencing this time from the time recorded on a receiver 
clock on GPS time. This time difference multiplied by the speed of signal propagation 
gives the pseudorange. Once code lock is attained, the code may be removed from the 
incoming signal and the satellite navigation message extracted. The signal, now demod- 
ulated, may at this point be used for processing carrier phase measurements. 

Squaring channel receivers operate by squaring the incoming satellite 
signal effectively removing any phase reversals resulting from code modulations im- 
pressed on the signal at the satellite. This produces a signal that is double the frequency 
of the incoming signal (See Figure 4). The phase difference between this squared signal 
and an internally generated replica is still easily determined however. The accumulated 
phase difference (whole and fractional cycles) is continuously counted and from this the 
carrier based change in range is determined by multiplying by the wavelength of the 
squared signal. 
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Squaring Channel 





Figure 4. Squaring Channel 


GPS receivers are equipped with quartz crystal oscillators rather than 
the more precise, stable, and expensive atomic oscillators used aboard the satellites. 
These oscillators do however provide more than sufficient precision and short term sta- 
bility for the typical length of a surveying session. Like the satellite clocks, the receiver 
clocks control signal generation. 

Internal receiver software estimates the receiver position by employ- 
ing least squares or Kalman filtering algorithms using the broadcast ephemeris satellite 
positions and pseudorange or carrier phase measurements. All receivers estimate the 
receiver clock bias in addition to the postion parameters and then use this value to reset 


the receiver clock to GPS time. 
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b. GPS Operating Principles - Emphasis on Point Positioning 
(1) GPS Observables and Observation Equations. There are two types of 
observation that may be exploited to determine positions, the pseudorange measurement 
and the carrier phase measurement. For point positioning particularly with GASP, the 
carrier phase is of primary interest and so 1s emphasized here. Pseudorange measure- 
ments to the extent they are used by GASP, serve only to monitor the validity of the 
carrier phase measurements and are not used in the actual position estimation. The 
fundamentals of pseudorange positioning are briefly touched on, however, in an effort 
to convey some general concepts and introduce some common terms. 
(2) Pseudorange. Under ideal conditions (1.e., precisely known satellite 
location, signal propagation speed, and perfectly synchronized satellite and receiver 


clocks) the travel time of the signal would be given by 
t=1,—1, 


where ¢, is the time of signal reception at the receiver and f, is the time of satellite signal 
transmission. This time difference would then be converted to distance units (p) by 
multiplying by the appropriate speed of signal propagation (c), where p=ct. This 
range measurement, along with the range measurements of two additional satellites and 
the known locations of the three satellites, could then produce a system of equations 
that can be solved to uniquely determine the receiver positions. 

In the real world, we must be concerned with error sources and their 
effects. A major source of error in any GPS measurement is the inability of the receiver 
and satellite clocks to maintain alignment with the reference time standard (GPS time). 
Of these two types of clock error, that associated with the receiver dominates due to the 
lower precision and stability of these clocks. The Allan variance is often cited as an in- 
dication of oscillator stability. The Allan variance values for receiver crystal quartz 
oscillators are typically several orders of magnitude larger than the satellite atomic 
oscillators. Satellite clocks take longer to warm-up and to stabilize initially but they 
maintain stability much longer than do receiver clocks. 

To address the problem of clock errors, we introduce a combined 
satellite and receiver, first order clock correction term as a fourth parameter in our sys- 


tem of equations. Expressing the receiver clock offset from GPS time as 


T,={,—1 
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where, 
t is the receiver clock time 
t is GPS time 


and the satellite clock offset from GPS time as 
Ta = li! 


where, 


t, is the satellite clock time 


the combined clock correction term becomes 
dT=T,—T, 


Because of clock and other types of errors, the measured range between satellite and 
receiver is not the true range but is rather a biased range referred to as the pseudorange. 


Designating the pseudorange as p, the system of observation equations is represented 
by: 


Pr, = Pit c di; 


with i greater than or equal to 4, 


where p, the true range from the ith satellite, is given by 


Pi =F [(xs, - sal ar Os, -= y) T (Zs, dE Ar E |X; = X,| 


In this equation, 
x, is the x coordinate of the ith satellite 
y,, 1s the y coordinate of the ith satellite 
z, is the z coordinate of the ith satellite and 


X, Yn z, are the unknown receiver coordinates. 


The addition of this fourth parameter requires a minimum of four observation equations 


to solve for the receiver coordinates and the clock correction. If t, is now introduced to 
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express the time delay associated with other error sources the pseudorange equation 
becomes: 


Pr, = Pi (aa Tp) 


The pseudorange is the observable most often used to obtain a navi- 
gation solution where submeter level accuracy is not required. Three dimensional navi- 
gation solutions currently offer best case positional accuracy at 16 meters. While 
certainly suitable for navigating a vessel, a position in error by 16 meters could not be 
used as a reference for establishing or extending mapping control. For high precision 
relative or point positioning, phase observables are needed. Because the carrier phase 
wavelength is shorter than the code modulations to the carrier signal (20 centimeters for 
Ll as opposed to 30 meters for the P-code), the carrier signal may be used to provide a 
more precise distance measurement. 

(3) Carrier Phase. It is possible to obtain distance (or more exactly 
change of distance) information by measuring the phase of the carrier signal. The dif- 
ference between the phase of a receiver generated carrier signal and the incoming, 


Doppler shifted, satellite carrier signal gives the carrier beat phase observable 


O meas = ais p + fat =$) 


where, œ is in cycles and fis the signal frequency. 

The total continuous carrier beat phase measurement consists of an 
accumulating count of the whole and fractional difference in cycles since the time of 
signal acquisition. What is not recorded is an unknown number of integer cycles at initial 
signal reception (the whole number of cycles between the receiver and satellite). The 
unknown number of cycles is referred to as the cycle or integer ambiguity. The total 


continuous carrier phase can be expressed 
P totat = Pmeas + N(%) 


where N(i,) is the integer ambiguity. A carrier phase equation comparable to the 


pseudorange equation can be written, 


Py, = À Pmeas, 
= p,— c (dti + Tp) + A N(to)i 
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where, p, is the phase biased pseudorange, 
t, is the signal time delay associated with other error sources and 


A=c/fis the signal wavelength. 


See Figure 5 for a graphic representation of the carrier phase concept. 
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Figure 5. Components of Carrier Phase Measurements 


c. Measurement Errors and Error Models 
There are two types of errors that effect GPS measurements and degrade the 
accuracy of the estimated position, systematic errors and random errors. Systematic 
errors are errors resulting from a predictable source. These errors are typically of con- 
stant magnitude under a set of given conditions. When the cause is understood, these 


errors can be removed by a correction model or a change in observation procedure. 


15 


Measurement biases resulung from systematic errors produce offsets between the values 
that are observed and the “true” value. After systematic errors are eliminated from 
measurements, random errors remain. These errors are defined by repeated measure- 
ments. A repeated measurement will not yield identical values each time due to the oc- 
currence of random errors. They indicate how well a measurement may be repeated and 
characterize the precision of a measurement. If they are the only error source, they are 
a measure of accuracy. Random errors are usually small and the probability of a positive 
or negative error of a given magnitude are the same. Random errors are dealt with in a 
least squares adjustment. [Ref. 7, 8] 

GPS systematic errors can be divided into three general categories: satellite 
errors, station errors, and observation dependent errors. Satellite errors include errors in 
the computed satellite ephemeris or in the model for the satellite clock provided in the 
navigation message. Station errors occur as a result of errors in the receiver clock. Ob- 
servation dependent errors are related to factors that influence the speed of signal 
propagation, ambiguities in phase observables, etc. 

(1) Satellite and Receiver Clock Errors. Inability to perfectly synchro- 
nize the receiver clock and the more tightly controlled satellite clock to a standard ref- 
erence (or GPS) time produces time offsets. If, for instance, there exists a one 
microsecond timing misalignment between satellite time and receiver time, a 300 meter 
range bias will result. Additionally, the clocks may be misaligned in frequency (a.fre- 
quency bias) or the frequencies may change over time (frequency drift). Any of these 
wil produce clock errors. Random error (or noise) also effects the time measurements. 

The timing behavior of the satellite clocks is carefully monitored by 
the ground stations and the drift of the clocks from standard GPS time determined. The 
amount of the satellite time offset, frequency offset, and frequency drift are then ex- 
pressed as coefficients of a second-order polynomial that are transmitted in the naviga- 
tion message. The polynomial coefficients are determined via a least squares adjustment 
performed at the master control station in Colorado Springs. The polynomial has the 
form, 


2 
Tx = 09 +4) (1— to) + a (t — to) 


where, &% is some reference epoch, 
æ is the satellite clock time offset, 


a, 1s the frequency offset, and 


16 





a, is the frequency drift. In this manner, satellite clock synchronization to GPS 
time is typically maintained to within 20 nanoseconds. 

GPS receiver clock errors must be modeled if the user desires a high 
accuracy solution. Receiver clock errors may be treated in a fashion similiar to that for 
the satellite clocks. The polynomial coefficients may be estimated as additional param- 
eters in a least squares adjustment along with the receiver coordinates. 

(2) Orbit Errors. The positions of the satellites with time are well known 
but, due to forces acting on the satellites that may not be adequately modeled, are not 
perfectly determined. The error in satellite position propagates to contribute to error in 
the receiver position. The magnitude of this error depends on whether the broadcast 
(predicted) ephemeris or the precise (post-fit) ephemeris is used to determine the satellite 
positions. 

Improving the models used to describe the forces acting on the satel- 
lite is one approach to obtaining more accurate satellite positions. Another is to include 
parametric models for the forces as part of the orbit estimation process performed by the 
control segment. Depending on the application, the effects of orbit errors may be ig- 
nored altogether or the data may be differenced to reduce or eliminate these effects. 

(3) Observation Dependent Errors. The speed of signal propagation is 
influenced by many factors including; signal interaction with the ionosphere and 
troposphere, relativistic effects on the signal, etc. 

(a) Signal Interaction with the Ionosphere— Signal interaction 
with the free electrons found in the ionosphere produces a change in path length that 
may be on the order of tens of meters. The free electrons are released from gas molecules 
ionized by incoming solar ultra-violet radiation. Any condition that acts to release more 
electrons, such as increasing the amount of incoming solar radiation (midday or during 
increased sunspot activity), will correspondingly lengthen the signal path. Satellite- 
receiver geometry also plays a role as the total number of free electrons along the path 
is a function of the distance that the signal travels through the ionosphere. Thus when 
a satellite is near the horizon the signal encounters more electrons than when near the 
zenith. This ionospheric effect is frequency dependent and is inversely proportional to 
the square of the signal frequency. A comparison of measurements on LI and L2 may 
be used to derive a dual frequency correction given by 


Los, ~~ P12] 
dpi; = E AA 
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where F is defined by 


| = ( Ju ) = | 
Si2 
For single frequency receivers, an ionospheric model must be used in place of the dual 
frequency correction or ignored altogether. [Ref. 2: p. 489] 
(b) Tropospheric Interaction— Tropospheric effects are a result 
of refraction in the neutral atmosphere and are not frequency dependent. This effect may 
be separated into two components: the dry component, and the wet component. The dry 


component comprises about 90% of the total effect and is a function of the surface at- 


mospheric pressure. It is approximated in the vertical by 
DTC = 2.27 x 10 “(m/mb) P, 


where P, (the atmospheric pressure) is in millibars and DTC (the dry term range con- 
tribution) is in meters. In the zenith direction, this corresponds to a maximum range 
bias of about 2.5 meters. The DTC increases with decreasing satellite elevation angle and 
at five degrees above the horizon it ranges from 20-30 meters. 

Estimating the effects of the wet component is a more compli- 
cated proposition. It depends on the total water vapor content along the signal path and 
hence on the temperature, pressure, and humidity. Surface temperature, pressure and 
humidity values are used to estimate the magnitude of the integrated effect by employing 
an atmospheric model such as the Hopfield or Chao models [Ref. 2: p. 490]. If obser- 
vations from water vapor radiometers are available, they will most accurately profile 
atmospheric conditions along the path. 

(c) Integer Ambiguity— In differential positioning applications, 
missing whole cycles (refered to as integer ambiguities) in the carrier phase measure- 
ments must be resolved if we are to fully exploit the more precise nature of phase 
measurements. Efforts must be made to account for any whole cycles not recorded in 
the observations. This may be the result of a loss of phase lock between receiver and 
satellite due to some obstruction to the signal. Also the initial number of whole cycles 
between satellite and receiver (the integer ambiguity) must be determined. Many tech- 
niques have been developed to resolve the absense of any cycles from the measurements. 
In point positioning with GASP, resolution of the integer ambiguity is not a problem 


as we shall see. 
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Other effects that must be accounted for include: relativistic 
effects, earth rotation (1.e., receiver position in motion) during time of signal trans- 
mussion, and offset of the satellite transmitting antenna position from the position of the 
satellite center of mass (the position given by the ephemeris). Corrections to account for 
the effects of these error sources are applied to the phase biased pseudorange and the 
pseudorange observables before the GASP model is formed and the estimation algorithm 
1s executed. Details on the models employed to compute these corrections are given in 
the chapter on GASP processing. 

(4) Comments on Selective Availability. Selective Availability (SA) is the 
intentional degradation of the position solution available to a select segment of the user 
community (1.e., most civilians and unfriendly military) operating in a real-time naviga- 
tion mode. It may be implemented by the DoD at their discretion. The solution degra- 
dation is accomplished by broadcasting inaccurate positions for the orbiting satellites 
or by dithering the satellite clock so that an inaccurate signal transmission time is ob- 
tained. This will not greatly effect relative positioning applications since receiver differ- 
encing schemes remove these errors. For point positioning applications that rely on the 
post-computed precise ephemeris, only the clock dither will effect the solutions. 

d. Geometric Effects 

In addition to the errors just presented, the configuration of the satellites 
relative to the receiver also has an effect on the determination of the receiver position. 
This geometric influence is referred to as the Dilution of Precision or DOP factor. It is 
roughly represented by the ratio of the positioning accuracy (o,) to the measurement 
accuracy (c) or 


DOP - 

Ge 
Actually there are a number of DOP factors. The most commonly referred to of these, 
the Geometric Dilution of Precision (GDOP), is more precisely defined as the square 
root of the trace of the covariance matrix (the trace being the sum of the four diagonal 
elements; the three position variances and the time variance). The DOP is a measure 
of the geometric strength of the satellite configuration that changes with time. The time 
of most favorable observation for obtaining a user position is when the DOP is small 
(with a DOP of less than five preferred) and falling. Since the number and the positions 
of the satellites visible to the observer changes over time, many possible geometric sce- 
narios are available over a tracking session. A GPS survey should be planned to take 
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best advantage of changing geometric scenarios. For point positioning, a geodetic 
quality position solution can usually be obtained with approximately four hours of data 
collection. During times of very favorable geometry, a geodetic quality point position 
may be obtained in a shorter time. 
e. Differencing Techniques 

By forming linear combinations of the basic pseudorange or carrier phase 
equations, a number of error sources common to the measurements being differenced 
will cancel or be greatly reduced. The accuracy of the computed position can be signif- 
icantly improved by employing differencing because many sources of measurement error 
are removed or reduced. The three types of differencing combinations often used are: 
differencing between two satellites, differencing between two receivers, and differencing 
between two epochs (time periods). These are referred to as single differences. The dif- 
ferencing schemes involving the carrier phase equations are the only ones presented here 
since they are used to form the GASP model. Differencing schemes that involve the 
pseudorange equation are ignored since pseudorange observations do not contribute to 
the position estimation except as a means of screening the carrier phase observations. 
Also, differencing techniques involving multiple receivers are neglected in this treatment 
except to state what types exist and the what errors they effect. 

Bétween-epoch single differencing involves differencing two equations of the 
instantaneous carrier phase for one satellite and one receiver. The carrier phase 
equation at epoch | is described by 


Py, = Pi — ¢ (at, + t)) + 4 M(t) 
The carrier phase equation at epoch 2 is given by 
Py, = Pa — ¢ (dt, + ty) + A No) 
Differencing then yields 
Pg, — Pg, = (Pi — P2) — ¢ (dt, — dt, + Tt) — tp.) 
Introducing 6 as the difference notation the equation may be written 


Opy = dp — c (ddr — 5t,) 


See Figure 6 for a graphic depiction of between-epoch single differencing. 
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Between—Epoch Single Differencing 


Satellite Q 


Figure 6. Between-Epoch Single Differences: Differencing Range Equations of 


One Satellite Over Two Consecutive Measurement Epochs 


The advantage of this scheme is that the initial integer ambiguity is removed. 
Similiar manipulations can be performed to yield between-satellite single 


differences (See Figure 7). For two satellites represented by the superscripts «a and $, 


phase equations are given by, 
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Po =p —c{dr + 18) + À Mio) 
and, 
ph = pf — c (drf +1) + À Nr)” 
The between-satellite single difference equation for the carrier phase is 
po — Pp =P" — p" — c (di — de + tp — Tp) +2 (M0) — No) 
or, 


dp, = dp + c (dtx — dtp) + A dN 


Between—Satellite Single Differencing 


Satellite 8 
Satellite © 


Figure 7. Between-Satellite Single Differences: Differencing Range Equations of 


Two Satellites Over One Measurement Epoch 
i | 
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For between-satellite single differencing, receiver clock errors are removed or reduced. 

It is also possible to form between-receiver single difference equations to 
remove satellite clock and orbit errors. This is important in relative positioning but 
cannot be used in point positioning. 

We can now form double difference equations from the single difference 
equations. For example, a satellite-receiver double difference may be formed by differ- 
encing two between-satellite single differences (involving the same pair of satellites) cver 
two receivers. This scheme removes or reduces the effects of errors associated with sat- 
ellite and receiver clocks. 

Receiver-time double differences and satellite-time double differences may 
be formed in a similiar manner. Receiver-time double differencing removes or reduces the 
effects of errors associated with satellite clocks and eliminates integer ambiguities. 
Satellite-time double differencing (essentially the GASP model, see Chap. III for more 
details) removes or reduces the effects of errors associated the receiver clock and elimi- 
nates integer ambiguities. 

It is also possible to construct an equation for a receiver-satellite-time triple 
difference observable. This gives the change in the receiver-satellite double difference 
from one epoch to the following epoch. With triple differencing, in addition to the can- 
cellation of integer ambiguities for carrier phase measurements, all clock and satellite 
orbit errors are removed. Triple differencing is widely used in relative positioning appli- 
cations [Ref. 9]. This technique cannot be employed in point positioning. 

The disadvantages to differencing are that the number of observations has 
been reduced and mathematical correlations are introduced as a product of the differ- 
encing process. This produces a weaker solution than that provided by not differencing. 
Correlation matrices should be computed to equate the differenced data to the undiffer- 
enced data but generally differenced observables are treated as uncorrelated. 

Now that systematic errors have been reduced or removed by error model- 
ing (1e., applying measurement corrections) or by differencing and the mathematical 
model has been defined by forming the undifferenced or differenced observation 
equation, the receiver coordinates and other parameters may be estimated by invoking 
the method of least squares. For a more complete explanation of the use of the least 


squares method in GPS, see Appendix B. 
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2. Positioning with Transit Doppler 

The method of Transit Doppler positioning is briefly reviewed here. Since this 
method was used to provide the independent reference positions to which the GPS de- 
rived point positions were compared, some insight into Transit Doppler positioning is 
necessary. 

The U.S. Navy Navigation Satellite System or Transit has been in continuous 
Operation since the mid 1960's. Transit can be regarded as the forerunner of GPS in 
many respects. Many of the ideas and techniques developed during the era of Transit 
positioning have been refined and employed in GPS. It is anticipated that Transit will 
be replaced by GPS in the near future. Like GPS, Transit consists of three segments: the 
satellites, the ground control, and the users. 

Six active Transit satellites are deployed in circular polar orbit at an altitude of 
approximately 1100 kilometers. This height was selected because excessive orbital height 
would provide too low a rate of change of Doppler frequency. If the orbit were too low, 
the Doppler frequency rate of change, ionospheric refraction, and especially the effects 
of atmospheric drag become too great. The orbital period of the satellites is about 107 
minutes. See Figure 8 for a general representation of the Transit Positioning System. 
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1. Receives 
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Transit Satellite Positioning 


Figure 8. The Transit System: Overview of the Transit Doppler System Seg- 
ments and Functions 
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Each satellite transmits two separate frequencies at 400 MHz and 150 MHz. 
The frequency generation is controlled by a single highly stable crystal oscillator. The 
use of two frequencies permits the determination of an ionospheric refraction correction. 
A broadcast message containing orbital information is superimposed on the two carrier 
frequencies by phase modulation. 

Three major differences between GPS and Transit may already be noted; the 
differences in satellite height, the use of only six satellites as opposed to 24 GPS satel- 
lites, and the use of crystal oscillators rather than the more precise atomic oscillators 
on-board the satellites. By having only six satellites available with fewer than this actu- 
ally visible, it takes much more time to acquire an equivalent number of Transit obser- 
vations. 

Tracking stations record the Doppler measurements (i.e., the Doppler shift in 
the frequency transmitted by the satellite) on each satellite pass. This information is re- 
layed to the central processing or control station where the satellite orbits are determined 
and extrapolated, then this ephemeris data 1s updated and uploaded to the satellites 
(about every 12 hours) for subsequent rebroadcast to the user. In addition, a timing 
station is responsible for monitoring the time signals received from the satellites and 
adjusting the satellite clocks as needed. 

a. The Basic Principle of Doppler Positioning 
The received frequency will differ from the transmitted frequency due the 
Doppler effect because the receiver and satellite are moving relative to each other. If the 


satellite transmits a stable frequency f, then the frequency at the receiver is given by 


f= 


where, 

r= A is the range rate 

r is the distance or range between receiver and satellite 

c is the speed of signal propagation. 

The Doppler frequencies (or Doppler shifts) are measured by subtracting 

the received shifted frequencies f, from a constant receiver reference frequency f, (the 400 
and 150 MHz frequencies alluded to earlier). The time of closest approach of the satel- 
lite is the time when f equals f. Integrated Doppler measurement techniques that count 


the number of accumulated cycles of Doppler shift are utilized in most Transit Doppler 
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receivers. The reason for this is that counting cycles can be performed more precisely 
than instantaneously measuring the frequency. 

The receiver position may be determined by continuously counting the 
number of cycles of the Doppler frequency. Referred to as the Continuously Integrated 
Doppler (CID) measurement mode, these observations along with accurate positions for 
the satellite can be used to establish the receiver position in two dimensions. For three- 
dimensional positioning, multiple satellite passes must be observed. 

The Doppler counts must be corrected for the effects of atmospheric re- 
fraction in both the ionosphere and the troposphere. In this respect, all other errors 
(relativity, earth rotation, etc.) that effected the GPS measurements will be present in the 
Transit Doppler measurements as well. Attempts, such as error modeling, should be 
made to minimize their impact. 

For the point positioning mode, the one of interest in this study, the ob- 
servations from multiple satellite passes were collected with a single Doppler receiver 
over about a four day period. Precise ephemerides were used to compute the satellite 
positions from which the CT coordinates of the site were determined. The estimated ac- 
curacies for a position solution achievable under this type of scenario is believed to be 
on the order of one meter. 
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IL EXPERIMENTAL PROCEDURE 


A. GENERAL 

The Monterey Bay Precision Positioning Experiment was conducted in early De- 
cember 1990. Aspects of the experiment relevant to point positioning will be emphasized 
in this and following sections. A large data set of static and kinematic GPS measure- 
ments was collected over four consecutive nights, from December 4 through December 
7, with four commercially available GPS receivers. Each night’s collection session lasted 
about six hours, from 11:00 p.m to 5:00 a.m. local time. This time window was selected 
because it was the period of maximum satellite visibility for the week of the experiment. 
This would permit the tracking of at least four satellites simultaneously, a crucial con- 
sideration from the standpoint of kinematic operations (at least four satellites are needed 
to solve for a three- dimensional position and clock bias). Fortunately, this was also the 
time when ionospheric effects on the measurements were minimal. All GPS measure- 
ments were recorded at a rate of one per second. 

The MBPPE static GPS positioning data was collected at the Naval Postgraduate 
School’s Beach Lab in Monterey, California. An array of five marks, refered to as the 
Doppler Array, was established to serve as the receiver reference points. In order to in- 
itialize the kinematic operations, a reference array (the Lobos Array) was established at 
the Monterey Bay Aquarium Research Institute’s pier in Moss Landing, California. One 
of the points in this array would provide additional static GPS measurements that would 
be utilized in this thesis. Another site, established on the roof of the NPS Mapping, 
Charting, and Geodesy building (in Monterey), would supply static GPS measurements 
for the Selective Availability test. See Figure 9 for a map of the experiment area. 
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Figure 9. Map of Experiment Area 


B. EQUIPMENT 
l. Acquisition and Familiarization 


The four models of GPS receiver used in the experiment were the TI 4100, the 
Trimble 4000ST, the Ashtech LD XII, and the Magnavox MX 4200. Two rubidium 
oscillators accompanied the TI 4100 receivers. For a summary of some of the features 


available on each, see the following tables (the Magnavox is excluded since it was not 


used in point positioning). 


M 1. ASHTECH —— oe FEATURES 








Standard Features | Features Comments 


No of Chanae 12 dual [12 dual frequency chan- chan- 


Internal RAM Capacity 6 ie 


Recording Interval 
Dual Frequency L1/L2 


L1/L2 Carrier Phase, C/A 









Table 2. TRIMBLE RECEIVER FEATURES 


Standard Features 
AA A RI AAN 
No. of Channels 8 dual frequency channels 
. 1 MB - used external 
Internal RAM Capacity datalogger (PC) 
Recording Interval 
Dual Frequency L1/L2 


Measurements L1/L2 Carrier Phase, C/A 
code Pseudorange on Ll 








Two of each type of receiver was acquired either from the vendors themselves 


or from some other source. The Ashtech, Trimble, and Magnavox models were bor- 


rowed from the respective vendors, while the TI 4100s were on loan from the Texas 


Department of Highways and Public Transportation. The receivers were obtained a few 


weeks prior to the data collection target dates so that the experiment participants could 
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Table 3. TI 4100 RECEIVER FEATURES 


Standard Features Comments 
E 
No. of Channels 4 dual frequency channels 


Internal RAM Capacity cae via 


cording Interval 
Dual Frequency LIL 


L1/L2 Carrier Phase, P 
Measurements code Pseudoranges on 
L1/L2 


become properly familiarized with the equipment. One person was assigned the respon- 
















sibility of learning the operational aspects of a particular receiver (as applicable to ex- 
periment requirements). Once acquainted with the equipment, a document describing the 
essential details of the receiver operation was drafted by each individual. Then all other 
participants were cross-trained on the different receivers. This would ensure that every- 
one would be able to start-up operations and trouble shoot if the situation arose. 
2. Equipment Set-Up and Data Collection 

One of each type of receiver was located at a static shore site, the Naval Post- 
graduate School’s Beach Lab. The others were located aboard the Research Vessel 
POINT SUR where the kinematic segment of the experiment was conducted. 

a. Beach Lab 

Before the commencement of data collection, the locations of the stations 

over which the individual receiver antennas would be set-up had to be established. At 
the Beach Lab shore site, the absolute position of a pre-existing mark (identified as 
DOP) had previously been established by the Transit Doppler method. This was used 
as the reference mark from which the positions of five new marks would be established 
(designated DOP1 - DOPS). This configuration of closely spaced shore marks (about five 
meters separation between marks) was refered to as the Doppler Array. The monuments 
were set and the positions of the marks were determined by employing conventional 
terrestrial survey techniques using a steel measuring tape to to obtain distances and a 


Wild T2000 theodolite to observe horizontal and vertical angles. 
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After the conclusion of the experiment data collection stage, the Transit 
Doppler station position for the mark DOP3 was resurveyed using a MX 1502 Transit 
Doppler receiver furnished by the DMA. Because the position for the mark had been 
established prior to the October 1989 Loma Prieta earthquake, a post-earthquake posi- 
tion for the mark was required. The Transit Doppler data was collected over four days 
and was submitted to DMA for the determination of the DOP3 position solution. The 
positions of DOP3 and the other stations comprising the shore site array were then up- 
dated based on the new information produced by this Transit Doppler survey. These 
were the independent position solutions to which the GPS solutions were compared. 

Each receiver-antenna was assigned to an individual monument. The 
Ashtech antenna was set up over DOPI, the Trimble over DOP2, the Magnavox over 
DOP4, and the TI4100 over DOPS. Except for the second night of data collection, when 
two antenna cables were inadvertantly attached to the wrong receivers (affecting the 
Ashtech and TI receivers), this configuration remained intact throughout the exper- 


iment. See Figure 10 for the Doppler Array layout. 


N 


N 


HORIZONTAL DIRECTION 
TAPED DISTANCE 
GPS LINE 


DOPPLER SCALE 1I cm = 2m  (1:200) 





Figure 10. The Doppler Array: The NPS Beach Lab Receiver Monument Sites, 
Ashtech Set-up Over DOP1, Trimble Over DOP2, TI 4100 Over DOPS 
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The Krupp Atlas Polartrack range and azimuth laser positioning system was set up over 
mark DOP3. This was used to determine the reference trajectory of the ship in the 
kinematic operations. 

The heights of the receiver antennas from their marks were measured both 
at the start and end of nightly operations. These measurements are necessary to adjust 
the solution from the antenna electrical center to the mark. 

Meteorological data (i.e., temperature, pressure, and relative humidity) was 
not recorded on site but was obtained from the NPS Meteorology Department for the 
nights of the experiment. Because the location of the School’s meteorology recording 
station is in close proximity to the Beach Lab site (less than 1000 meters), the weather 
data should closely reflect conditions at the Beach Lab. 

b. The Lobos Site 

In order to initialize the kinematic operations, the pre-departure and post- 
arrival positions of the POINT SUR, GPS antenna array (located atop the mast in the 
crows nest) had to be established. The positions of four reference points located in the 
Vicinity of the POINT SUR's dock were determined for this purpose. This station group 
is refered to as the Lobos array (See Figure 11). 
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Figure 11. The Lobos Array: The Moss Landing Monument Sites, TI 4100 Re- 
ceiver Set-up Over LOBOS3 


Following the experiment, the position of one of the points, LOBOS3, was 
resurveyed with both GPS and Transit Doppler. The TI 4100 receiver was used for the 
GPS survey. The GPS data was collected continously over a three day period (February 
5-7, 1991) at a measurement collection rate of 30 seconds. This supplied three inde- 
pendent data sets, one for each day of the survey. The Transit Doppler survey, con- 
ducted from January 31 - February 4, 1991, established the Transit Doppler position for 
the mark. The independent GPS solutions are compared to the Transit Doppler solution 


in later sections. 
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The height of the TI 4100 antenna above the mark was measured at at the 
beginning and at the completion of the three day collection session. Weather data was 
again obtained from the NPS Meterology Station. Although the degree to which this 
data accurately represents conditions at Moss Landing is not known, it was felt it should 
adequately reflect LOBOS3 weather conditions. 

In determining the point positions for the Beach Lab and Moss Landing 
sites, the precise satellite ephemerides and clock models were obtained from DMA for 
the time of applicability. This provided a comparison of the solutions estimated with this 
information versus the solutions estimated with the broadcast ephemeris and clock data 
collected during the surveys. 

c. Building 224 

As of early July 1991, Selective Availability (SA) had been activated by the 
DoD. GPS data was collected for two days (July 16 and 17) at a site that had been 
previouly established on the roof of the NPS Mapping, Charting, and Geodesy center 
(site: BLDG 224.3). The data was collected with two different TI 4100 receivers at a 
measurement rate of 30 seconds. The antenna height above the mark was recorded. No 
weather data was obtained for this test. The data was then processed with GASP using 
the broadcast ephemeris only and the site position solution was compared to the known 
site position in order to evalute the effect of SA on the computed position. 
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IHI. DATA PROCESSING WITH GASP 


A. PROCESSING FLOW OVERVIEW 

Before the raw receiver data could be processed, some preliminary operations had 
to be performed to convert the data into a format suitable to the GASP program. The 
raw data was translated on a PC into a GASP compatible ASCII format (i.e., FICA or 
RINEX). The actual procedures and programs involved depend upon the particular re- 
ceiver and are covered in more detail in the next section. Once the data had been con- 
verted and all other required data assembled in the proper formats, the processing was 
initiated. 

GASP consists of two major program units. The first, called the GeoSTAR 
PREProcessor or STARPREP, accepts input from a FICA or RINEX data file and from 
files containing satellite ephemeris and clock information, meteorological data, and sta- 
tion information (such as antenna height above mark and a priori station coordinates). 
STARPREP then computes and applies a series of measurement corrections to the L1 
phase biased pseudorange and pseudorange measurements. Once corrected, these meas- 
urements are referred to as two-frequency corrected measurements. Also, measurement 
time tags are adjusted to reflect the correct GPS time of signal transmission and the 
satellite positions are interpolated to correspond to the corrected time tags. This infor- 
mation, along with the corrected observations and the values computed for the meas- 
urement corrections, is output to a “point” file. A station information file is also 
generated containing the a priori station position, antenna height, and other relevant 
data. These files are subsequently used as the input to the second major program unit, 
GASP. 

The GASP unit accepts input from the point and station files, forms the GASP ob- 
servable, and performs the estimation for the point position in CT coordinates. See 
Figure 12 for a general depiction of the processing flow. 
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Figure 12. GASP Processing Flow 


B. PERSONAL COMPUTER (PC) OPERATIONS 
The operations required to convert the receiver data into a GASP acceptable format 
varied for each receiver. An outline of the general procedure for each receiver is pre- 
sented. 
l. TI 4100 Receiver Data Conversion 
The raw tracking data collected with the TI4100 Basic External Processor Pro- 
gram (or BEPP) operating system had to be converted to the Floating point, Integer, 
Character, ASCII (FICA) format. This involved the use of utility programs developed 
by the University of Texas, Applied Research Laboratory (UT/ARL). The raw tracking 
data was converted to the Floating point, Integer, Character (FIC) binary format by the 
program GS2FIC. This format was then converted to the FICA format by the FICFICA 
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program. Because the Beach Lab data was collected at a one second measurement rate, 
it Was necessary to decimate the data to 30 seconds in order to alleviate file storage 
problems. Also, GASP documentation suggests a 30 second measurement rate and over 
our collection periods it would provide more than a sufficient number of observations 
to produce meter level position solutions. 
2. Trimble Receiver Data Conversion 
The first step in the conversion of the Trimble raw tracking data to the RINEX 
format was to download the data collected by the receiver into a restructured binary 
format. TRIMVEC (Trimble supplied processing software) accomplished this task. Then 
the TRRINEX programs, developed by the creators of the RINEX format (a team from 
the University of Bern, Switzerland), were utilized to convert the Trimble binary format 
to the RINEX format. The Beach Lab, Trimble data was also decimated to 30 seconds. 
One difference between the FICA and RINEX formats should be noted at this point, the 
broadcast ephemeris and clock information obtained by the receivers from the Naviga- 
tion Message is presented as separate data blocks in a single FICA data file. In the 
RINEX format, this information is presented as a separate file altogether. This fact be- 
came important When considering the RINEX modifications to GASP. 
3. Ashtech Receiver Data Conversion 
The Ashtech raw data also had to be downloaded before RINEX conversion 
could be implemented. The Ashtech GPPS processing software performed both the 
downloading and the RINEX conversion. ASHTORIN was the program used to per- 
form the conversion. It too permitted the Ashtech Beach Lab measurements to be 
decimated to 30 seconds. 
4. PC to VAX Data Transfer 
The observation and broadcast ephemeris data, now in either the FICA or 
RINEX format, had to be transferred to the NPS Digital VAX computer. This is where 
the GASP programs and program code resided and where processing would be done. 
The data was transferred to the VAX via 9-track tape using the Overland Data Tape 
Software, DEPOT program. The precise ephemerides and clock information, furnished 
by the DMA for the weeks of applicability, were also transfered to the VAX. 
Meteorological and station data files were created in the GASP specified formats for the 
different days and receiver stations. See Figures 13 and 14 for examples of the 
meteorological and station files. 
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352800. 
356400. 
360000. 
363600. 
367200. 
370800. 
374400. 
378000. 
381600. 
385200. 
388800. 
392400. 


DANDNNAROO ru 
220000000090 


Field 1 - 2-digit year identification 
l'icld 2 - GPS week number 
J-icld 3 - GPS time of weck 


licld 4 - Barometric pressure in millibars 


Field 5 - Temperature in Celsius 
Field 6 - Relative humidity in percent 





Figure 13. Sample STARPREP Meteorological File Contents and Record De- 


scription 
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Record | A 

Field 1 - Site 1.D. (forinat 18) 

Field 2 - Record type (format Al) 

Field 3 - Station 1.D. (format F10.0) 

Field 4 - Station type (format Al) 

l'icld 5 - Station name (forinat ^30) 

Field 6- Station Latitude in rads (format D16G.10) 
Field 7 - Station Longitude im rads (format D16.10) 
Field 8 - Station Latitude in degrees (format D15.9) 
Field 9 - Station Longitude in degrees (format D15.9) 


Field 10 - Station height in meters (format F10.6) 
Record 2 sá 

Vicld 1 - Site I.D. 2 (format 18) 

Field 2 - Record type (format Al) 

Field 3 - Site type (format Al) 

Field 4- Antenna offset north of mark (format 18.6) 
Field 5- Antenna offset cast of mark (format F8.6) 


Field 6- Antenna height offset (forinat 8.6) 
Ficld 7 - Year, day of antenna set up (format 15) 
Ficld 8 - Elevation angle cutoff for 16 sv's (format 161"4.1) 


Figure 14, Sample STARPREP Station File Contents and Record 


C: 


Description: Sv’s, referred to in record 2, field 8, is another term for 


satellites. 


DATA FILE EDITS AND RUNSTREAM GENERATION 
Prior to STARPREP processing, the following information contained in the FICA 


or RINEX observation files was validated: 


the a priori station coordinates 
the antenna height 
the year, month, and day of session 


week-crossovers | E 


e the PRN numbers of the satellites tracked 


e the measurement record interval (or rate) 


Special attention had to be paid to the last item since the decimation process did not 
automatically update the measurement rate in the header block for the TI 4100 data. 
Much of the data collected at the LOBOS3 site, corresponding to periods when 
fewer than three satellites were available, was discarded. This still provided three inde- 
pendent data sets (one for each day of site occupation) of between eight and 11 hours 
duration at a 30 second recording rate. The BLDG 224.3 data set, collected contin- 
uously over two days at 30 seconds, was divided into five independent subsets ranging 
from about four to nine hours. The data sets and subsets, time spans, and the number 


of satellites used from the individual data sets are identified in the following tables. 


Table 4. BEACH LAB DATA SET IDENTIFICATION AND 
TIME SPANS 


: Time Span No. of Satellites 
Tigo | 62 | 8 O 


338 | Ashtech | Se |" 9 | 
A Ll o O 


Ti4i00 | 62 _ | 8 O 
339 | Ashtech | 6 | _9 
IC E a 


AA AA EA 
340 | Ashtech | 68 | 9 | 
Timbe | 75 | 9 


TI 4100 | 65 | 8 
341 | Ashtech | 69 | 9 
_ Trimble | 67 | 9 O 
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Table 5. LOBOS3 DATA SET IDENTIFICATION 
AND TIME SPANS 


Time Span (hours) No. LES 
een ee ee TS | 


36 1 A 
a ools || 8 
n= ee 










Table 6. BLDG 224.3 DATA SET IDENTIFICATION AND TIME 
SPANS 


Da Time Span No. of Satellites 
y (hours) Used 





A runstream on the VAX was generated that signified what data files were to be 
used as input, output or temporarily created, the unit numbers assigned to each file, 
which of the available satellites were to be used as the source of observations, whether 
broadcast or precise ephemeris and clock states were to be used, the corrections to be 
applied, etc. The runstream is a VAX batch file that begins execution of the STARPREP 
routine, deletes temporary files at the completion of processing and, in general, controls 
the processing operations. See Figure 15 for a sample precise ephemeris runstream. 
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Figure 15. 


command procedure to run STARFREP (PRECISE ETREMERIS ANO CLOCKS) 
The GPS data flle (FICA format) aust be loaded In 

the directory: (MALYS.STARFILES}. The Precise Ephemerides 

and Precisa Clock Files must elso ba located In [MALTS.STARFILES) 
VRITE SISSOUTPUT ” 1 AM RUNNING STARPREP. STANDBY, IT WON'T TAKE LONG.” 
set default [bredthauer.star{lles) 
define eyeSoutput ster346ltPl.out 
dir 

type [bredthaver.starprep.OBJECT)runstar341RP,con 
$ run [bredthever.eterprep. OBJECT Jatarprep3 
ASGFILES 04 STATIONLOG. OAT 
ASGFILES 07 RINDOP2343t .0AT 
ASGFILES 10 PTDr2t0341P.PNT 
ASGFILES 11 DOFI341A.90N 
ASGFILES 12 O2BLBAOJ41T.TEN 
ASGFILES 13 O6BI.BAO343T. TEM 
ASCFILES 14 O9BLBAOQIAIT.TEM 
ASGFIJ.ES 15 33DI.BAO0I43T.TEH 
ASGFILES 16 1281.8A0341T. TEN 
ASGFILES 17 1585.8A0341T. TEN 
ASGFII-ES 18 1681.8A0341T. TEN 
ASGFILES 19 18B1.BAQI4IT.TEM 
ASGFILES 20 19BLBAOI41T. TEM 
ASCF3IES 24 EF290336.EPR 
ASGFILES 25 EF690336.EPR 
ASGFILES 26 EF990336.2PA 
ASGFILES 27 €F1190336.Et8 
ASGFILES 26 8F1290336.8PA 
ASGFILES 29 EF1590336.RFR 
ASGFILES 30 EF3690336.2th 
ASGFIIZS 31 EF1090336.EPR 
ASGFILES 32 EF1990336.RPR 
ASGFILES 36 MTdop20341P.0AT 
ASGFILES 37 RCdop2034)P.O0AT 
ASGFIIES 386 OVdop20341P.DAT 
ASGFILES 39 STOP2t0343P.DAT 
ASGFILZS 40 7090336. EFh 
ASGFILES 41 MMET361.DAT 
ASGFILES 42 MRCV.DAT 
ASGFILZS 43 MSAT.DAT 
ASGFILES 44 MSTA341.0AT 
ASGFILES 31 TEMPO2.TEM 
ASCFILES 932 TPHPO6.TEN 
ASGFILES 353 TENFD9.TEN 
ASGFILES 34 TENPF)).TEN 
ASGFILES 53 TENP)2.TEN 
ASGFIIES 36 TEMPIS.TEM 
ASGFILES 57 TEMPIG.TEM 
ASGFILES 50 TEMP10.TEM 
ASGFILES 59 Tenri9. TEN à 
PTROCSEO EDIT TTACCOR OATACOR 
CHPUTCDR ID IR ER TH SC CR SA 
ATPLTCOR 10 IR ER TH SC CR SA 
PLOTCORR 10 IR ER TH SC CR SA 
EFHEMRIS P 
EDCNTROL * RCVOPT TOLOPT METOPT 
DFOUCGSON  OBUG 
LOGICALS SETERK 


SELECTSY 2 6 9 11 12 15 16 18 19 
ENOINFUT 


OELETE 4,TEM/* 
l a o o o e e o e DR SES a a e e o en ceccceesescoesons § END OUT UT PILES TO CASP 


$ rename/log pt0r2t0341P.pnt {bredtheuer.gasptiles]* 
$ rename/log st0r2t0341P.det [bredthever.gaspilles)]* 
$ desos syeSoutput 
$ 
$ 
$ 


CA LP LR LR APR LA LP LR LA 2 
> == = == 
9 


-=THAT'S ALL THERE IS TO IT 





i 
vrite eys$output “STARPRE? ÑAS COMPLETED. See STARO20.0UT for results.” 


Sample STARPREP Runstream for Precise Ephemeris 
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D. STARPREP 
I. Input 

Much of the information presented in the following sections on STARPREP 
and GASP processing is found in Malys, et al [Ref. 1]. 

STARPREP accepts input data from an observation file (pseudorange and car- 
rier phase measurements in either FICA or RINEX format), broadcast or precise 
ephemeris and clock data from file input (or in the case of the FICA broadcast 
ephemeris data as part of the FICA observation file), a meteorological file, and a station 
file. The carrier beat phase measurements are converted to kilometers by multiplying 
by the nominal L1 or L2 wavelength. In the TI 4100 FICA observation file, the time of 
signal transit between satellite and receiver is given rather than the actual pseudorange. 
This is converted to the pseudorange by multiplying by the speed of light in a vacuum. 
RINEX presents this as the pseudorange originally so that no conversion to distance 
units 1s necessary. 

2. Error Models 

The time and data corrections applied during the course of STARPREP proc- 
essing will be briefly discussed and the error models used to compute the corrections 
presented. This material is originally presented by Malys, et al [Ref. 2: pages 489-491]. 

a. Time of Transmission 
Measurement time tags are adjusted from time of reception (¢,) to time of 


transmission (f,) by, 


where, p, 1s the pseudorange observation at t, and 
c, the speed of light in a vacuum, equals 299792.458 kilometers per second. See 
Appendix C for a list of the constants used in the GASP programs. 
b. Satellite Clock 
For the broadcast ephemeris, the predicted satellite clock parameters avail- 
able in the navigation message are used to compute the satellite clock offset fram GPS 
time by, 


2 
Tx = @ + (tx — lo) + a, (1, — to) 


where, 4, is the time of applicability of &, a,, a 
æ is the predicted satellite clock time offset 


| 


a, is the predicted frequency offset 
a, is the predicted frequency drift. 


For the precise ephemeris, the precise clock states are used to compute the satellite clock 
correction at the observation epochs. 
c. Receiver Frequency Offset 

The TI 4100 receiver LI and L2 reference signals are offset from the LI and 
L2 carrier frequencies by -6000 Hz and + 7600 Hz respectively. A correction is computed 
which removes the number of cycles in the data due to these biases. The epoch of the 
initial recorded carrier beat phase measurement for a satellite is used as the reference for 
subsequent carrier phase measurements. For example, the correction for L1 phase data 


in cycles is given by, 
Nec = — 6000 (1 — t,) 


where, N,,. 1s the number of cycle counts 
tis some affected epoch, in seconds, and 


t is the initial epoch recorded for the satellite, in seconds. 


The Trimble and Ashtech receivers do not have offsets in their reference signals so the 
values assigned to the offsets were set equal to zero for RINEX processing runs with 
these receivers. 
d. Ionosphere 
Signal interaction with the free electrons found in the ionosphere produces 
a change in signal path length. This ionospheric refraction index is frequency dependent 
so a comparison of measurements on LI and L2 may be used to derive a dual frequency 


correction. Defining the quantity F as 


[4-1 


a dual frequency ionospheric correction for the pseudorange is given by, 


(Pr, P Pr) 
aD, > F 
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This can only be applied to data collected with the TI 4100 receiver since it is the only 
model that supplies pseudoranges from the L2 signal. 

The dual frequency ionospheric correction for the carrier phase for the TI 
4100 receiver is given by, 


(Po. 000 
dp, = RTS + Lau (— 6000) — 412 (7600)] (A‘) 


where, À, and 1, are the transmitted wavelengths and Ar is the interval between the initial 
epoch of observation and the epoch being corrected [Ref. 2 : p. 489]. For the Trimble 
and Ashtech receivers with no receiver frequency offset, this reduces to, 


A (Op aoe) 
Po, — le 


e. Troposphere 

Tropospheric effects are a result of refraction in the neutral atmosphere and 
are not frequency dependent. This effect may be separated into two components: the dry 
component, and the wet component. The dry component comprises about 90% of the 
total effect and is a function of the surface atmospheric pressure and satellite elevation 
angle. Estimation of the wet component is more difficult than estimating the dry com- 
ponent. It depends on the total water vapor content along the signal path and hence 
on the temperature, pressure, and humidity. 

In GASP, one of two tropospheric models may be selected to compute the 
integrated tropospheric correction, the Chao or Hopfield models. Both models use sur- 
face weather data and the satellite’s elevation angles to compute the correction. The 
temperature, pressure, and relative humidity are used to compute the wet and dry com- 
ponents of the zenith tropospheric delay values (Z,,,, Z,,,). The satellite elevation angles 
are then used to compute the wet and dry multipliers (F,,, F.,, ). The total correction, 
for either model, may be generally given by, 


dP rop TF (Zary Fay Te Lwet iter 


See Chao and Hopfield (Refs. 10, 11] for more detailed explanations of these two models. 

The Hopfield model was selected as the tropospheric model for all the re- 
sults produced in this study. A few of the similiarities and differences between the two 
models should be noted. Both models produce similiar zenith values for the wet and dry 
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components. The major difference is the dependency of the correction on the satellite 
elevation angles. The Chao model, originally developed for use in arid locations, is not 
as accurate as the Hopfield model in estimating the wet tropospheric component off the 
Zenith [Ref. 12]. 
f. General Relativity 
If using the broadcast ephemeris, the general relativity correction is com- 


puted as a function of the broadcast orbital elements and 1s given by, 
MP 
do, =20 | ME Jio a'h E sinn] 


where, GM, the product of the universal gravitational constant and the Earth's mass, 
equals 3.986005 x 10%km?*] sec?, and 

e is the eccentricity of the satellite orbit 

ais the semi-major axis of the satellite orbit 


E, is the eccentric anomaly at the observation epoch. 


For a more thorough description of the broadcast orbital elements see Appendix A. 
If the precise ephemeris is used, the correction is given by, 
DO 

AP re = = 

where, X is the satellite position vector and 
V is the satellite velocity vector at the observation epoch. See [Ref. 13] for a full 
explanation of this effect. 
g. Earth Rotation 


To account for the fact that the Earth is rotating while the signal is traveling 


from satellite to receiver, the Earth rotation correction is given by, 


dpe = + LX — Xn) Xi — (Xi — Xn) X% 


where, w, the WGS 84 value for the Earth’s rotation rate, equals 7.2921151467 x 10-5 
radians per second, 
X is the interpolated satellite position vector 


X, is the receiver’s a priori position vector 
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1, 2 are the vector components along the X and Y axes of the WGS 84 reference 
frame. 
h. Satellite Antenna Offset 
The precise ephemeris contains the positions of a satellites center of mass. 
To adjust this position to the electrical center of the satellites transmitting antenna, this 


satellite antenna offset is approximated as, 
dpsa=R°S 


where, R is the range vector from the a priori receiver position to the satellite and 
S is the scaled vector from the satellite's center of mass to the satellite sub-point. 


S is given by, 
S = SAo (8) 


where, SA, is the L-band satellite antenna offset in the nadir direction (0.88m) and 
e, is the unit vector from the center of mass of the Earth to the satellite. 

The broadcast ephemeris gives the position for the phase center of the sat- 
ellite transmitting antenna [Ref. 14: p. 83]. This would make application of this adjust- 
ment unnecessary When using the broadcast ephemeris. The GASP program does not 
distinguish between the precise and broadcast ephemeris in this regard. This results in a 
range error that is common to all satellites when the broadcast ephemeris is used. Be- 
cause this error is common to all satellites, the GASP algorithm will effectively remove 
this error. 

The data sets were originally processed with this correction applied for both 
ephemerides. Two of the previously processed data sets Were reprocessed using the 
broadcast ephemeris, this time without the satellite antenna correction applied. The 
maximum difference for any component over the two data sets was on the order of two 
centimeters. It can be reasonably concluded that the GASP differencing scheme has 
eliminated this source of error. 

3. Output 
The STARPREP preprocessor generates three main files; the point file, the sta- 
tion file, and the output file. GASP accepts the point and station files as input. The 
output file contains plots of all the applied data corrections for each satellite over time, 
plots of the uncorrected range and carrier phase measurements over time, and a sum- 
mary of file information. This file is helpful in identifying problems that may have arisen 
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during the course of processing. The point file contains the two-frequency corrected 
observations, the values computed for the individual corrections, the corrected meas- 
urement time tags, and the interpolated satellite positions (in CT coordinates) corre- 
sponding to the corrected time tag. The station file contains the a priori station 
coordinates and the height of the antenna above the mark. See Figure 16 for a general 


summary of STARPREP processing. 


Station end Satellite Ephemerides 
Meteorological -_—_ 5 end Clock States 
info 


Corroute and Apply 
Corrections 


Adjust Time Tags 


Dual Frequency Corrected 
A Priori Station 
Position end Deta. Tive Taga 


Anterra Height and e 





Figure 16. STARPREP Processing 


E. GASP 
1. GASP Runstream Generation 

Prior to GASP processing, a GASP runstream was created. This runstream 
specifies the station and point files to be used as input, assigns a filename to the GASP 
output file, and if desired, permits the selection of non-standard GASP processing 
options. Standard options preset certain processing values or features but these may be 
overridden by the operator. See Table 7 for a list of the available processing options and 
standard values and features. 
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Table 7. GASP RUNSTREAM PROCESSING OPTIONS AND 
STANDARD VALUES OR FEATURES 


Processing Options Values or 
Features 
[Plotto Printer À PP 


Standard Deviation on A Priori Position Components 

0.05 
(km) 
PRN Number as ‘Base’ Sat, Default is 00 for Sequenc- RUES 
Ing 


Number of Batch Least Squares Iteration 
Number of Sequential Estimation Iterations — 

















2. Automatic Data Editing 

After the point and station files have been entered into GASP, two kinds of data 
editing are performed to ensure consistency between the pseudoranges and phase biased 
pseudoranges observed from each satellite. First, the corrected pseudoranges over two 
successive epochs are subtracted and the corrected phase biased pseudoranges corre- 
sponding to the same epochs are also subtracted. If these “delta” ranges differ by more 
than some user specified tolerance, the GASP observables (i.e., the two carrier based 
phase biased pseudoranges) are rejected from the data set. This is referred to as 
pseudorange editing. Note that this is the only use of the corrected pseudoranges. They 
serve only to monitor the phase biased pseudoranges for outliers and are not used in the 
estimation process. 

For the processing runs performed for this study, the standard pseudorange ed- 
iting tolerance of five meters was initially selected. After a few runs with the Trimble 
and Ashtech receiver data, it became apparent that a significant percentage of the data 
was being rejected from the data set (on the order of 35% or more). In order to avoid 
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this, the pseudorange editing tolerance was relaxed to ten meters. This held down the 
amount of data rejected to under 20 percent in most cases. There are two possible ex- 
planations for why such a large amount of data was rejected. The first, and most prob- 
able, an ionospheric correction could not be computed for the Trimble and Ashtech 
pseudorange data since no pseudorange measurements were supplied on the L2 signal 
for these receivers. This would produce some discrepancy between the corrected 
pseudorange and the phase biased pseudorange where an ionospheric correction was 
applied. A second possible explanation is that the pseudoranges are observed from the 
C/A-code modulations rather than the more precise P-code modulations. The increased 
measurement noise on the C/A-code observations may also have contributed to the 
higher rejection figures. 

The second edit test performed on the data uses the Root-Mean-Square (RMS) 
of the residuals of the previous iteration as the rejection criterion. The RMS is initialized 
before the first iteration. Any phase biased pseudorange observables with a residual 
greater than three times the RMS of the previous iteration are rejected from the data set. 
These editing schemes are executed prior to the formation of the GASP observable. 

3. The GASP Observable 

The between-epoch single difference equation is used to form what is known as 
the GASP observable. Two consecutive carrier based phase biased ranges from the same 
satellite are differenced. Recall that this between-epoch single difference equation .was 
given by, 


Ôp = Ôp —c (ddt — 6t,) 


This between-epoch difference is then differenced with the corresponding between-epoch 


difference from another satellite. The resulting equation is given by 


Aôp, = ôpy — ôph 
= (60° — 6p}) — c (57% — 87% + 679 — $14) 


where a and £$ signifies the two satellites. It may be expressed more conveniently as, 


Ap, = Adp — c (Adt, — Adt,) 


Since the satellite clock errors and the errors associated with various other sources have 


been modeled, the equation becomes, 
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Adp = Adpy + c (Adt, —AÔr,) 


The value Adp is the GASP observable. See Figure 17 for a conceptual representation 
of the GASP observable. 


The GASP Observable 


Satellite 8 
Satellite © 





Figure 17. The GASP Observable: Differencing Over Satellites Over Two Con- 
secutive Measurement Epochs 
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This type of differencing scheme reduces the effects of errors in satellite clocks 
and orbits, and removes receiver clock errors and integer ambiguities. For a four chan- 
nel receiver with all channels simultaneously tracking satellites, three GASP observables 
can be formed for a given epoch pair. At least two satellites per epoch pair are required 
to form an observable for that epoch pair. For each epoch pair, one satellite is used as 
the reference from which the others are differenced. The reference satellite is selected 
sequentially, so that for every new epoch pair processed, the next higher satellite PRN 
number is used as the reference. The selection sequence cycles back to the lowest PRN 
number once the list of tracked satellites is exhausted. Individual epochs are used only 
once to form an observable. 

Reference satellite sequencing reduces correlation among the observables. Pre- 
cision, the RMS of the residuals, and the variance-covariance matrices of the estimated 
parameters are all improved using this type of satellite selection rather than selecting the 
satellite with the most stable clock as the reference. [Ref. 1: p. 24] 

For example, if the Ll phase biased pseudoranges are represented by p, , four 
satellites simultaneously tracked signified by superscripts «, $, y, and y (in order of as- 
cending PRN number), and four consecutive epochs represented by subscripts 1, 2, 3, 


and 4, a between-epoch difference can be formed for each satellite given by 
Ph, — Pha = Para 
Ph, — Pha = $0, 
Ph, — Pa, = Po 
Po, — Pb. = Por, 


The three GASP observables formed by satellite sequencing for this first epoch pair are 
given by, 


a f aß 
Pr - OP $1, - Adps,, 

a Ve ee, ay 
OP¢,, T ÔP h E Aôp,, 


a n _ an 
OP 4, ™ Pb, 7 Adpg,, 
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Then using $ as the reference satellite for the next epoch pair at 3 and 4, the GASP 


observables are given by, 


B a pa 
Phi 7 Pr = Adpg,, 
B y = By 
Pr En OP; i A0p5,, 
B y Bn 
OP F Ph me Aóp5,, 


The next epoch pair at 5 and 6 would use satellite y as the reference and epoch pair 7 
and 8 would use satellite y. Reference selection would then return to satellite « for epoch 
pair 9 and 10. 
4. Batch Least Squares and Sequential Estimation 

After the GASP observables have been formed, a least squares technique is used 
to estimate the station position components in the CT coordinate system. These are the 
only parameters estimated. A clock correction parameter is not estimated. It is assumed 
that the modeled satellite clock states have sufficiently aligned the satellite clocks to GPS 
time and that the receiver clock error has been removed by differencing. The interpo- 
lated satellite positions are held fixed in the estimation providing the reference frame in 
which estimation takes place. After three iterations of a batch least squares, the esti- 
mated parameters, the variance-covariance matrix for the parameters, and the RMS of 
the residuals provide the input to the second estimation step (for more on the least 
squares method see Appendix B). 

This 1s a sequential estimation algorithm based on a Kalman filtering routine. 
The Kalman state is the receiver position vector. It does not include any clock states. 
Fundamentally, a Kalman filter updates measurements from one observation epoch to 
the next. This allows the parameter estimates and covariances to be updated at each 
measurement epoch. The batch method provides a measure of data noise to the se- 
quential processor. The final RMS of the residuals from the batch processor is taken as 
the variance of one GASP observable processed through the sequential processor. Since 
the sequential processor updates the station coordinates for each new observable proc- 
essed, the estimated station coordinates can be plotted as a function of time. Plots of 
the covariances and convergence in CT coordinates as a function of time may also sup- 
plied. 

For the processing runs performed in this study, a comparison was made be- 
tween the final batch least squares component estimates and the final component esti- 





mates that had passed through the sequential processor. The magnitude of the mean 
component differences, averaged for all processing runs using the precise ephemeris, was 
about 15 centimeters. 
5. Output 

The most important information contained in the output file generated by 
GASP is the estimated station position. This is presented as the X, Y, and Z coordinates 
of the CT coordinate system and in geodetic coordinates; latitude, longitude, and height 
(D, A, h) relative to the WGS-84 ellipsoid. Estimates for the uncertainty in the coordi- 
nates is obtained from the variance-covariance matrix for the estimated parameters. A 
correlation matrix is computed that expresses the linear independence among the esti- 
mated parameters from the parameter variance-covariance matrix. A posteriori standard 
deviations for the estimated station coordinates are obtained by taking the square root 
of the diagonal elements (i.e., the variances) of the variance-covariance matrix. These 
provide estimates for the precision of the estimated point position components. The 
GASP processing sequence 1s graphically depicted in Figure 18. 
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Figure 18. GASP Processing 


F. PROGRAM MODIFICATIONS 

In order for GASP to accept the data collected with the Trimble and Ashtech re- 
ceivers, it was necessary to modify the program to allow the introduction of data pre- 
sented in the RINEX format.2 This meant altering some existing program units and 
creating two new program modules for STARPREP. The first reads the RINEX obser- 
vation file. The STARPREP subroutines BLK006 and BLK124 that read the FICA 
measurement and station data blocks were used as models in building the RINEX sub- 
routine. The second module, called RINNAV, reads the RINEX navigation file that 
contains the broadcast ephemeris information. It was modeled after the STARPREP 


2 RINEX data may be presented in one of two possible forms; the original RINEX format 
and the RINEX 2 (version 2) format. GASP processing runs utilized data presented in the original 
RINEX format. Program code designed to process data in the RINEX 2 format, although present, 

bi is untested. [Ref. 15] , 
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subroutine BLK009 that reads the FICA data blocks containing the broadcast ephemeris 
data. 

To differentiate between FICA format and RINEX format input, a file naming 
convention was implemented. The program path to read either format is keyed to the 
first four characters of the input observation filename. For the RINEX data collected 
with the Ashtech receiver, the first four characters of the filename must be RINA. For 
the RINEX Trimble data, the first three characters must be RIN and the fourth char- 
acter anything other than A. For FICA data collected with the TI 4100, the first three 
characters may be anything except RIN. The reason for distinguishing between the 
RINEX Ashtech and Trimble data is due to a problem that occurred during processing. 
A wavelength factor of one was specified for the L1 signal and a wavelength factor of 
two specified for the L2 signal in the RINEX observation file for both the Ashtech and 
Trimble receivers. Using these wavelength factors with the Ashtech receiver caused all 
the observables to be rejected by the pseudorange edit. This was because of an incom- 
patibility between the internal receiver software and the receiver hardware. Apparently, 
the receiver software had not been updated to reflect changes in the hardware. To cir- 
cumvent this problem, the L2 wavelength factor for the Ashtech receiver is reset to one 
based on the RINEX filename. Once the receiver software is updated, the program 
statement that resets the wavelength factor will need to be deleted. 

Also keyed to the filenames are values assigned to the receiver frequency biases. The 
TI 4100 receiver incorporates receiver frequency biases of -6000 Hz on the L1 signal and 
7600 Hz on the L2 signal. The Ashtech and Trimble receivers do not incorporate these 
biases. 

Because the Trimble and Ashtech receivers possess more tracking channels than the 
TI 4100 and could thus supply simultaneous observations from more satellites, we felt 
we should take advantage of this feature by modifying GASP to accomodate data col- 
lected on these additional channels. The Trimble receiver records observations on eight 
dual frequency channels and the Ashtech receiver on twelve. Because increasing the di- 
mensions of the program arrays might cause space problems and adversely affect proc- 
essing, a cautious approach was taken and storage expanded to accomodate eight 
channels per epoch. Basically, this was accomplished by simply redimensioning the ar- 
rays that stored the input measurements, two-frequency corrected observations, time 
tags, and interpolated satellite positions. 

The program revisions to include more channels led to modifications to expand the 
maximum number of satellites that could supply observations over the course of a col- 
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lection session. The previous limit of eight was changed to twelve. Thus, the number 
of satellites allowed per collection session should not be a limiting factor to the potential 
benefits of using more receiver channels. If twelve satellites per session Were to be uti- 
lized, allowances for additional files associated with the extra satellites (and the RINEX 
broadcast ephemeris file) would be necessary. The maximum number of permissible 
program files was reset from 60 to 62 and program statements assigning unit numbers 
to particular files or setting limits on the unit numbers for certain file groups (such as 
the precise ephemeris and temporary files) were respecified. Specifically, the unit number 
assigned to the RINEX broadcast ephemeris file is 11, the upper and lower limits on the 
unit numbers for the first set of temporary satellite files changed to 12 and 23, limits on 
the precise ephmeris files changed to 24 and 35, limits on the output data files changed 
to 36 and 39, and the limits on the last group of temporary files changed to 51 and 62. 
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IV. RESULTS AND ANALYSIS 


In order to provide some measure for the accuracy of the computed positions, the 
known station position components are differenced from the GPS CT estimated position 
components for the various receivers and collection sessions. The known station posi- 
tions are the Transit Doppler positions. This supplies the AX, AY, and AZ offsets from 
the “true” positions. Component differences AX, AY,and AZ are converted to AE, AN, 
and AU (in local, cartesian east, north, and up coordinates) because this is a more fa- 
miliar and easily comprehended reference frame. The following discussion and analysis 
wil concentrate on the tables and plots that present the results in the local east, north, 
and up coordinate system. For the reader interested in the results presented in CT co- 
ordinates (X, Y, and Z), see the tables in Appendix D. 

The difference results are displayed in the form of tables and associated target plots. 
In the target plots, the Transit positions are represented by the origin and the symbols 
depict the values computed for the component differences. The GASP estimated formal 
error (one sigma) standard deviations are presented as the + terms in the difference ta- 
bles or as error bars in the difference plots. 

The known station positions are the Transit Doppler derived positions. Mark BLDG 
224.3 was established by a GPS relative positioning survey using the Transit Doppler 
mark DOP3 as the reference. See Appendix E for a summary of the estimated positions 
and estimated standard deviations on the positions for both the Transit Doppler and 
GPS point determinations. 

Repeatability, also defined as precision, yields a measure for the consistency of a 
given set of position results. The repeatability is determined by calculating the mean and 
the standard deviations (i.e., the observed errors) on a set of position component differ- 
ences. Two types of repeatability were determined. One averages over days by receiver 
accentuating the differences between the receivers. The other averages over receiver by 
day accentuating the difference between days. In the repeatability target plots, the mean 
differences are represented by the target symbols and the standard deviations are re- 


presented by the error bars. 


A. BEACH LAB RESULTS 
Table 8 presents the accuracy results of the Ashtech, Trimble, and TI 4100 receivers 
for the Beach Lab collection sessions using the broadcast ephemeris for the satellite po- 
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sitions. It can be seen from this table and the corresponding target plots (figures 19 and 
20) that the GASP estimated components give fairly good agreement with the Transit 
Doppler estimates. The average Root-Sum-Square (RSS) of the component differences 
(1.e., the magnitude of the difference vector) is at about the two meter level. The origin 
is within the ensemble of the position differences and formal errors. The estimates 
produced with the TI 4100 receiver were slightly better overall than the Trimble and 
Ashtech estimates but all were comparable. The Trimble results are the most widely 
dispersed and the Ashtech results slightly less scattered. The TI 4100 results exhibited 
the least dispersion of the three receivers. 

It also appears that, between the differences computed in the north and east di- 
rections, the greatest dispersion is produced in the east, 1.e., the estimated east compo- 
nents are less precise than the north. For all three components, the north component 
shows the least dispersion i.e., the best agreement is between the GASP estimated north 
component and the Transit Doppler north component. Comparison of the AH and AU 
component differences shows more dispersion in the vertical than in the horizontal. 

A lack of overlap between the GASP estimated formal errors, due to the scatter of 
the component differences, is also conspicuous. This suggests that the GASP formal er- 
rors underestimate the actual observed errors. The GASP formal errors appear to 
underestimate error in the east component to a greater degree than the north compo- 
nent. The dispersion of the vertical versus the horizontal differences suggests, overall, 
the GASP estimated errors were better at representing the horizontal errors than the 
vertical. This visual interpretation is supported by the values computed for the bottom 
line of the table which shows the mean component differences and standard deviations 
on the means over all receivers and days. The numbers indicate that the east and up 
components were underestimated by a factor of about three and the north by a factor 
of approximately two. 

Also, as compared to the differences computed for the other three days, the day 339 
results seem “off”. This may be due to poorer broadcast ephemeris satellite positions for 
this day. 


Table 8. BEACH LAB COMPONENT DIFFERENCES (AE, AN, AU) USING BROADCAST 
EPHEMERIS: Solution Differences Between GPS and Transit Doppler Methods, Sec- 
ond Term is GASP Estimated Formal Error, RSS is Magnitude of the Difference Vector 
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Figure 19. AEast versus ANorth for Beach Lab, Broadcast Epliemeris Data 
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Figure 20. AHor versus AUp for Beach Lab, Broadcast Ephemeris Data 
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Table 9 and Figures 21 and 22 present the accuracy results of the three receivers 
using the precise ephemeris. It can be seen that the position estimates produced with the 
precise ephemeris give better agreement to the Transit Doppler positions than the 
broadcast ephemeris. Overall, the precise ephemeris results are more accurate and pre- 
cise in every component difference as demonstrated by comparison of the mean com- 
ponent differences and the standard deviations computed about the mean differences. 
The overall RSS is at about the 1.5 meter level. Again, the origin is within the ensemble 
of all position differences and formal errors. 

The estimates produced by the three receivers are all comparable, with the Ashtech 
and Trimble giving somewhat better results than the TI 4100 receiver over the four days 
of observations especially in terms of difference dispersion. The component differences 
for the Ashtech receiver are all fairly tightly clustered about the origin while the Trimble 
and TI 4100 receivers show greater dispersion especially in the east component. The 
TI 4100 receiver also shows the greatest dispersion in the north and up directions. 

As seen in the broadcast ephemeris results, comparisons of dispersion between the 
north and the east reveal that dispersion is greater in the east than in the north direction. 
Also conspicuous is the greater dispersion in the vertical than in the horizontal. 

Again, we see that the GASP estimated formal errors have underestimated the ac- 
tual error as determined by the dispersion of the differences. For the precise ephemeris 
estimates, GASP formal error underestimates the actual error more in the east and up 
components. It represents the error in the north component fairly well. The actual error 
on the east and up components is underestimated by a factor of about two. The re- 
lationship for the error in the north component is almost one-to-one. 


Table 9. BEACH LAB COMPONENT DIFFERENCES (AE, AN, AU) USING PRECISE 
EPHEMERIS: Solution Differences Between GPS and Transit Doppler Methods, Sec- 
ond Term is GASP Estimated Formal Error, RSS is Magnitude of the Difference Vector 


Year 1990 AE (m) AN (ım) AU (m) AH (m) ‘| RSS (m) 


TI A100 ET 

338 “06+ FFE 
TI 4100 Dire 0357 
0747 


T1 4100 -16+.5 | —00+.5 | —1.6+.5 1.6 + .6 
340 —0.0+.6 | —04+.5 | -2.3+.5 FE 
—0.4 + .5 0.2+.4 | —1.3+.4 04+.6 | 14 | 































TI 4100 084-5 | 164.5 
341 024.6 | 064.5 | —01+.5 
1.4 + 6 0.5 4.5 


Mean Component Differ- 
ence and Standard Devi- —0.1 + 1.2 0.1 + .6 —0.7+ 1.1 1.1 +.7 1.6 + .7 
ation about Mean 
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Table 10 and Figures 23 and 24 display the Beach Lab repeatability results that av- 
erage over days by receiver for both the precise and broadcast ephemerides. The means 
and standard deviations on the means show that geodetic-quality point positions are 
achieved using the precise ephemeris. This is especially true for the results produced with 
the Ashtech and Trimble receivers. The standard deviations 1.e., observed errors are 
typically within or around one meter on all components for these receivers. The accura- 
cies, as determined by the RSS of the mean component differences, are within one meter 
for all receivers demonstrating that agreement to the Transit Doppler estimates are 
within the computed noise levels. 

The TI 4100 precise ephemeris results are less precise than either the Ashtech or 
Trimble precise ephemeris results and, interestingly, are less precise than the TI 4100 
broadcast ephemeris results. The reason for this is not known. However, examination 
of the RSS of the mean component differences reveals that the overall accuracy of the 
TI 4100 precise ephemeris estimates is better than the broadcast ephemeris estimates. 
This is presumably due to increased systematic error associated with the broadcast 
ephemeris satellite positions and clock states. 

The broadcast ephemeris estimates show that geodetic-quality point positions are 
not obtained using the broadcast satellite positions and clock states. Higher values 
computed for the RSS of the mean differences show reduced accuracy compared to the 
precise ephemeris estimates. Also, the broadcast ephemeris estimates are generally more 
imprecise than the precise ephemeris estimates. As previously mentioned, this was not 
true for the TI 4100 receiver however. 

It is also clear that averaging the component differences over days has significantly 
improved the accuracy and precision for all receivers for both ephemerides. The indi- 
vidual components show improvement and, as a consequence, the RSS of the averaged 
components is also improved. This demonstrates the importance of averaging over many 


independent position estimates to suppress the effects of random error. 
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The error bars, representing the actual observed error (the standard deviations about 
the mean component differences), show significant asymmetry. This dramatically illus- 
trates the reduced precision in the east component as compared to the north component 
and in the vertical component when compared to the horizontal. It also illustrates that 
the GASP estimated formal errors substantially underestimated the observed or true er- 
rors. For the precise ephemeris estimates, errors on the north component for all three 
receivers displayed close agreement to the GASP estimated errors. The Ashtech receiver 
showed the best agreement in the east component and the Trimble receiver in the up. 
The broacast ephemeris estimates also displayed close agreement in the north for all re- 
ceivers, but the east and up components showed even poorer agreement than the precise 


ephemeris. 
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Table 10. BEACH LAB REPEATABILITY (AE, AN, AU): AVERAGE OVER DAYS BY 
CEIVER: Mean RSS is Magnitude of the Mean Difference Vector 


en) 
(m) (n) (m) (m) (m) 


Precise 


T14100_ | 0.0+10 | -11+4 
0.72.1 1.9412 | 0.7+2.1 
0.8 +1,5 “19418 | 11+14 


TI 4100 0.3 + 1.7 —0.5+15| 0.6+1.5 
0.24.6 | —08+11] 054.7 
= 0.213 0.3 + 1,2 
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Figure 23. Mean AEast versus Mean ANorth for Beach Lab Repeatability: Average 
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Table 11 and Figures 25 and 26 display the Beach Lab repeatability results that av- 
erage over receivers by day for both the precise and broadcast ephemerides. The table 
and the plots show that the day 339 results, produced using the broadcast ephemeris, 
were the most inaccurate and imprecise of all the daily estimates. Overall, the broadcast 
ephemeris results showed much less precision than the precise ephemeris estimates. The 
precision on the positions estimated with the precise ephemeris was roughly two times 
better than those estimated with the broadcast ephemeris. This demonstrates the in- 
consistency of the broadcast ephemeris in estimating geodetic-quality point positions. 

Again, the plots show asymmetry in the magnitudes of the error bars. The east 
component displays greater imprecision than the north component. Comparison of the 
vertical and horizontal components does not reveal any obvious overall asymmetry ex- 
cept in a few individual cases. This also shows that, again, the GASP estimated formal 
errors underestimate the true errors. The computed values demonstrate, for the precise 
ephemeris, the north and up formal errors best estimate the observed errors for all days. 
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Table 11. 


BEACH LAB REPEATABILITY (AE, AN, AU): AVERAGE OVER RECEIVERS| 


DAY: Mean RSS is Magnitude of the Mean Difference Vector 


ones [De Mi an a | Maa 
P (m) (m) (m) (m) 


Broadcast 


Precise 


CEE Jon 14215] 0328 
Ss | 49220 | 0215 [19431] 19422 
orar | 004.7 | 10413] 0227 
3a conne] 014.6 | 0324 | 028.8 
[338 | -06zi7] 04410 | 1045 | 07419 
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B. LOBOS3 RESULTS 

The TI 4100 receiver was the only model used to collect data at the LOBOS3 site. 
Independent data sets were created for each of the three days of site occupation. Results 
produced by the broadcast and precise ephemeris are compared for each of the inde- 
pendent data sets. 

Table 12 and Figures 27 and 28 displaying the LOBOS3 component differences for 
both ephemerides clearly demonstrate the advantage of using the precise ephemeris to 
determine the point position. Accuracy is improved in every component. The broadcast 
ephemeris estimates show much lower levels of accuracy and precision and there appears 
to be a bias present in the south-east and up directions. Errors in the predicted satellite 
positions or in the broadcast satellite clock models are the probable cause of this bias. 

The precise ephemeris results show very good agreement with the Transit Doppler 
estimates and can definitely be considered geodetic-quality. As with the Beach Lab re- 
sults, the precise ephemeris estimates show greater dispersion in the east component 
than in the north. Any difference in dispersion between the vertical and the horizontal 
does not appear to be significant for either ephemeris. Also noteworthy, the TI 4100 
gave much better precise ephemeris results at this site than at the Beach Lab and the 
broadcast ephemeris results are much worse overall. 

Because of the peculiar TI 4100, precise ephemeris results obtained at the Beach 
Lab, it would be desirable to compare the precise ephemeris results obtained with the 
TI 4100 receiver at LOBOS3 to those obtained by the other two receivers at the Beach 
Lab. Direct comparison of the results requires some caution, however, because of the 
longer lengths of the data sets processed for LOBOS3. With this caveat in mind, com- 
parison of the LOBOS3 TI 4100 estimates to the Ashtec and Trimble Beach Lab esti- 
mates reveals that better results, in terms of both accuracy and precision, were produced 
with the TI 4100 receiver for all components at LOBOS3. This may be due to better 
precise ephemeris satellite positions provided for this particular week. Again, it should 
be stressed, although the LOBOS3 precise ephemeris results with the TI 4100 were best 
overall, all three receivers produced geodetic-quality point positions with the precise 
ephemeris. 
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Table 12. LOBOS3 COMPONENT DIFFERENCES (AE, AN, AU): Solution Differences! 
tween GPS and Transit Doppler Methods for Broadcast and Precise Ephemerides, 
ond Term is GASP Estimated Formal Error, RSS is Magnitude of the Difference Ve 


ee AE (m) AN (m) AU (m) AH (m) RSS ( 


CE Tire | -1625 | 1725 | 2128 | 21 
Broadcast [37 | 3448 | aste | 494.6 | 57410 | 73 
£ 


—0.4+.5 0.6 + .5 0.6 + .8 0.9 
Precise 0.7 
1.3.6 0.1 +.5 0.2 +.5 1.3.7 1.3 


78 


2 
dE PRECISE EPHEMERIS 
met 
z, 0 J Toe 
<j E 

DAY 
19 4 ee X 36 
° 37 
O 38 
=4 
-6 
22 0 


BROADCAST EPHEMERIS 


DAY 
X 36 
o 37 
O 38 





ER) 


Figure 27. AEast versus ANorth for LOBOS3 


ya! i 


79 


PRECISE EPHEMERIS 


AU (m) 


19} X 36 
o 37 
O 38 


BROADCAST EPHEMERIS 


AU (m) 


DAY 
X 36 
o 37 
O 38 





A 


Figure 23. AHor versus AUp for LOBOS3 


ly | D. nu 


80 


The LOBOS3 repeatability results averaging over days, displayed in Table 13 and 
Figures 29 and 30, clearly shows the inaccuracy and imprecision of the broadcast 
ephemeris results. The south-east-up bias alluded to earlier can also be clearly seen in 
the broadcast ephemeris estimates. Comparison of the error bars on the averaged precise 
ephemeris estimates shows the characteristic east versus north asymmetry. No signif- 
icant vertical versus horizontal asymmetry is evident, however. The magnitudes of the 
error bars on the precise ephemeris estimates show that, in this case, the GASP esti- 
mated formal errors more closely represented the true errors. However, the formal error 
estimates produced for the broadcast ephemeris results underestimated the true errors 
by a factor of approximately two for the east component, three for the north, and four 
in the vertical. The precise ephemeris statistics demonstrate that by averaging over in- 
dependent estimates, the accuracy and precision can be significantly improved. 


C. REGIONAL AGREEMENT BETWEEN GPS AND TRANSIT DOPPLER 

In order to establish a regional reference for the expected agreement between GPS 
and Transit Doppler position estimates in the Central California area, the overall accu- 
racy and precision of the position differences using the precise ephemeris are evaluated. 
The precise ephemeris mean component differences and standard deviations over all 
Beach Lab and LOBOS3 data sets were computed to provide values for this expected 
agreement. These values, obtained by averaging over 15 data sets, are mean AE 
= 0.21 + 1.1, mean AN = 0.09 + .52, and mean AU = —0.54 + 1.03. 


81 


Table 13. LOBOS3 REPEATABILITY (AE, AN, AU) AVERAGE OVER DAYS FOR TI 4 
RECEIVER: Mean RSS is Magnitude of the Mean Difference Vector 


Mean AË (m) | Mean AN (m) | Mean AU (m) | Mean AH (m) Me 


2.0 + 1.2 —3.2 + 1.5 4,1 + 2.2 3.8 + 1.9 5.6 + 2.9 
053 
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D. POSITIONING RESULTS - SELECTIVE AVAILABILITY ACTIVATED 

Table 14 and Figures 31 and 32 present the solution differences between the coor- 
dinates previously established for mark BLDG 224.3 and the coordinates estimated for 
the mark using the broadcast ephemeris after the activation of Selective Availability. 
From comparison of the broadcast ephemeris results obtained at the Beach Lab and 
LOBOS3 sites to the results displayed here, it is quite obvious that the activation of SA 
has severely impacted the position estimate. Both the accuracy and precision of the re- 
sults have been adversely affected. 

The average RSS of the estimated positions over the five data sets is about 16 me- 
ters. When contrasted to the approximately six meter average RSS computed for the 
LOBOS3 site using the broadcast ephemeris, some indication of the magnitude of this 
effect is realized. The east and up components are most affected. The north component 
seems to be much less effected by SA than the other two components. Precise 
ephemerides were not obtained for this survey period so a comparison of the broadcast 


ephemeris results to the precise ephemeris estimated position cannot be performed. 
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Table 14. BLDG 224.3 COMPONENT DIFFERENCES (AE, AN, AU): Selective Availabi 
Test - GPS Solution Differences Between Previously Determined Coordinates and | 
ordinates Determined After Activation of S/A (Using Broadcast Ephemeris), RSS 
Magnitude of the Difference Vector 


Day of 
ON AE (m) AN (m) AU (m) AH (m) RSS (n 


SA 23.64 1.1 | 20.5410 (2237005 31.3 
T 12.0+.8 | 3411 12.5 


C | 8820 | 2.7510 9.2 + 2.2 9.3 
198 | OD | -54+8 | 23+9 | 178+13 | 59+1, 18.7 
7.9 + 1.0 7.9 + 1.0 8.0 
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AHor versus AUp for BLDG 224.3, Broadcast Ephemeris Data 
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V. CONCLUSIONS 


GASP estimated point positions, produced with data collected from three GPS re- 
ceivers over seven days, are compared to independent Transit Doppler estimated posi- 
tions. The differences between the positions provide a measure of accuracy for the 
computed GPS positions. Both broadcast and precise ephemerides were used to estimate 
the GASP point positions. The precision for a group of estimated positions is deter- 
mined by calulating the mean of the differences and standard deviations on the means 


for the group. Major conclusions are: 


l. GASP point position estimates generated with the precise ephemeris agree to 
Transit Doppler point position estimates to within one meter on each axis. 


2. All three GPS receivers used in this study are capable of producing point positions 
with observed errors on the RSS of less than 1.5 meters, when the precise ephemeris 
and clock states are utilized. 


3. Point positions estimated with the precise ephemeris are more accurate and precise 
than positions estimated with the broadcast ephemeris. 


4. Selective Availability has had a dramatic effect on positions. 


The GPS point positions estimated by GASP using the precise ephemeris and clock 
states demonstrated good agreement with positions estimated by the Transit Doppler 
positioning system. The level of agreement was within the one meter noise level estab- 
lished for both systems. In general, the Transit Doppler positions were within the en- 
semble of corresponding GPS position estimates. This would seem to indicate that many 
of the systematic errors that effect the GPS measurements had been successfully modeled 
or removed by the GASP algorithm. Of course, systematic errors common to both po- 
sitioning systems may remain. 

Accuracy and precision indicate that all three GPS receivers examined in this study 
are capable of producing geodetic-quality point positions. Between-receiver compar- 
isons of the estimated positions using both the precise and broadcast ephemeris reveal 
that the Trimble and Ashtech receivers provide results that compare favorably to those 
produced by the TI 4100 receiver. This is not surprising considering that all three re- 
ceivers record the carrier phase measurements essential to GASP processing and high- 
precision point positioning. Although the carrier phase observable recorded by the three 


receivers is the same, the pseudorange observable exhibits some important differences. 
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Pseudoranges are observed from the LI signal using the C/A-code for the Ashtech 
and Trimble receivers. The TI 4100 receiver observes pseudoranges on both the Li and 
L2 signals using the P-code. Consequently, a dual frequency ionospheric correction can 
not be computed for pseudoranges observed with the Ashtech and Trimble receivers. 
This, along with higher noise levels associated with C/A code measurements, caused a 
high incidence of data rejection in preliminary processing runs. Relaxing the 
pseudorange edit tolerance to ten meters allowed the acceptance of a much higher per- 
centage of observations. Even so, rejection figures for these receivers were still higher 
than the TI 4100 receiver. 

Though a higher percentage of observations were rejected with the Ashtech and 
Trimble receivers, a high degree of accuracy and precision was maintained. The high 
accuracy and precision levels observed with these receivers may be attributed to their 
increased channel capacity. The ability to collect and subsequently to process more si- 
multaneous satellite carrier phase measurements provides greater geometric diversity. 
More measurements per epoch and the incorporation of more satellites into the proc- 
essing run provides increased system redundancy. These two factors act to strengthen 
the estimated position solution. Position solutions produced by the TI 4100 receiver 
should also be improved now that more satellites are allowed per processing run. The 
benefits gained from these additional measurements will be nullified if too many obser- 
vations are rejected in the data editing process. Thus, it may be desirable in future tests 
to determine an optimal pseudorange edit tolerance. 

Point positions estimated using precise ephemerides and clock states are superior to 
the solutions produced with the broadcast ephemeris and clock states. The results ob- 
tained with the three receivers at the Beach Lab demonstrate that, overall, better posi- 
tions estimates were provided with the precise ephemeris and clock states. This is clearly 
the case at the LOBOS3 site where significantly better results were achieved with the 
precise ephemeris. These results, together with results produced in past studies, demon- 
strate that positions produced with the precise ephemeris are generally more reliable, 
accurate, and consistent than those produced with the broadcast ephemeris [Ref. 2: p. 
497]. 

The recent reactivation of Selective Availabilty has accomplished its intended pur- 
pose. For the results reported in this study, SA has severely degraded the accuracy and 
repeatability of point position estimates that were produced with the broadcast 
ephemeris. The effects of SA on position estimates using the precise ephemeris were not 
assessed in this investigation. This may be a topic for future consideration. 
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The RINEX modifications to the GASP program now allow it to utilize data col- 
lected with a greater variety of GPS receivers. This will permit point position determi- 
nations from a multitude of additional sources and will enable DMA to augment its 
point positioning data base. Only data that was converted to the original RINEX format 
was used in these tests. Additional tests should be conducted with data presented in the 
RINEX 2 format to ensure that the program code designed to process this data is free 
of error. 

Additional tests should be conducted to verify the validity of the conclusions pre- 
sented here. Tests, similiar in design but extending over longer time periods (weeks as 
opposed to days), would provide a more satisfactory indication of the long-term repeat- 
ability for the different receivers. 
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APPENDIX A. KEPLERIAN ELEMENTS 


Six Keplerian elements and a reference time are necessary to completely describe the 
satellite orbit and the position of the satellite in the orbit. The six Keplerian elements 
given in the broadcast ephemeris are: 

e Right ascension of the ascending node (Q) 

e Inclination (1) 

e Argument of perigee (w) 

e Semi-major axis of the elliptical orbit (a) 

e Eccentricity of the orbit (e) 

e An element describing the position of the satellite on the orbital ellipse (e.g., Mean 

anomaly (M)) which is a function of time 

The ephemeris reference time is designated (z,,). The five Keplerian elements that de- 
scribe the orbit are Q, i, w, a, and e. M gives the position of the satellite in its orbit at 


a time t. See Figure 33 for a graphic representation. 
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Figure 33. The Keplerian Orbital Elements: Symbols are Explained in the Text 


Other Broadcast Ephemeris parameters that are transmitted in the Navigation Message 
describe the deviations of the satellite motion from the smooth ellipse defined by the six 
Keplerian elements refered to above. These perturbing terms include: 


e Mean motion difference (An) 
e Rate of right ascension (Q) 


T + Rate of inclination (/) i, F h 
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e Corrections to argument of latitude (C,, C,,) | 
e Corrections to orbital radius (C,, C,,) 
e Corrections to inclination (C, C) 
In order to transform the Keplerian orbital parameters to coordinates in the con- 
ventional terrestrial system, the following computations must be performed: 
e find time since reference epoch 
lk =I loe 


e solve for mean anomaly at t, 


My= My +| -= + An | 
a 


where the gravitational constant 
3 
u = 3.986005 x 1014 2 
sec 
e solve Kepler’s equation iteratively for eccentric anomaly, E, 
M, = Ex — sin Ex 


e compute true anomaly, 


i= an} (1 — e?) EA ue | 


(cos E, — e) 
e compute the argument of latitude, 

up = w + fg + Cue COS 2(w + fg) + Cus sin 2(w + fy) 
e compute the orbital radius, 

r, =a(l1—ecos E,)+ C,, cos 2(w + fy) + C,, sin 2(w + fy) 

e compute the orbital inclination, 

i, = lo + ty + Ci cos (w + fy) + Cy sin (w + fy) 
e compute the longitude of the ascending node, 

Ay = % + (Q — w.) ly — We loe 
where the mean earth rotation rate is 


w, = 7.292115147 x 10-* 224 


e then compute the CT (earth-fixed) coordinates, 
Xp = rp COS Up COS Àp — rg SÌN Uy COS i, Sin A, 
Yp = rp COS up Sin Àg + rg Sin up COS i, COS Ay 


Ly = Tr sin 4; sin lp 


94 


APPENDIX B. LEAST SQUARES IN GPS 


The contents of this appendix follow the development in Uotila [Ref. 16]. This de- 
velopment focuses on the general, nonlinear, observation equation model. This is the 


model utilized by GASP in performing its batch least squares adjustment. 


A. THE MATHEMATICAL MODEL 


Using matrix notation the mathematical model is given by, 


L* = F(X") 
eee 
L? — £= KX) 


where, Le = Theoretical values of the observed quantities 
X= Theoretical values of the paramters 
Le= Estimates or adjusted values of the observed quantities 
X= Estimates or adjusted values of the parameters 
Lt = Observed values of the observables 


€ = “True” errors 


B. OBSERVATION EQUATIONS 
Prior to employing the least squares method, a nonlinear function is linearized by 
expanding in a Taylor series and eliminating terms second order and higher. For a set 


of nonlinear functions in matrix form we get, 


L -e= FAX) = LX) + LE yen gt (X=) t = FU) + Se |e. Xt 


OF 
A rex 


where, X° are approximate values for the parameters and is the point about which ex- 
pansion is done, and 


A is called the design matrix and it has the form 
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OX Sf Oe 
df of, Of 
A == (¿DOS 
Ya M À 
ôx, ôx, ` ôx, 


where, n is the number of observations, and 
u is the number of equation parameters. We will assume that redundancy exists, 
i.e., the number of observations is greater than the number of unknown parameters. 


Using the relations: 
c=AX+FX°)- L! 
V AE ee 


where, V is a vector of differences between the observed values and the estimated values 
(V is refered to as the residual vector) and 


X=X -Xx 
KEK 
and from, 
F(X°) = L° 
Coan 


where, L° are the computed observations evaluated at X we get, 


A 
V=AX+L 
These are the observation equations. 


C. MINIMUM VARIANCE SOLUTION 

To find a unique set of parameter estimates some condition must be imposed on the 
residuals. By minimizing the sum of the squares of the weighted residuals, we obtain a 
best estimate for the parameters in a least squares sense. 
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The weight matrix denoted P is given as, 
P=0, EL. 


where, Z;! is the inverse of the observation variance-covariance matrix and 
o, is the a priori variance of unit weight. 

The weight matrix provides a way to assign relative importance to various observa- 
tions. The observation variance-covariance matrix lets us express random errors or 
noise associated with the observations as variances and relationships between observa- 
tions as covariances. It is through the variance-covariance matrix that knowledge about 
the observations is propagated. This allows the variances and covariances for the esti- 
mated parameters to be determined (i.e., it gives an estimate of the noise in the solution). 
The variances make up the diagonal elements of the matrix and the off-diagonal ele- 


ments are the covariances. The variance-covariance matrix may be written, 


2 
Or, O11 12- + > FL Ln 
2 
Ot2 11 %12 + + + FL2Ln 
2, = 
2 | 
On L1 OLn L2: e . O Ln 


where, o, = o}, i.e., a symmetric matrix. The variance-covariance matrix is diagonal if 
the observations are uncorrelated. 

For the minimum variance solution, we want to minimize V7 P V which is the sum 
of the squares of the weighted residuals in matrix notation. By some algebraic manipu- 


lation we can show, 
VTPV=XATPAX+2L PAX+LTPL 
Minimizing by partially differentiating with respect to X gives, 


1 êV" PY) 


> ——= 4 PATA P 
ax 


This is known as the normal equation. Now solving for X we get, 
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Ÿ=—(4 PA)! ATPL 


where, X is the adjustment vector that is applied to X° giving an improved estimate for 
the parameters. A number of iterations may be necessary before X reaches some ac- 


ceptable level, that is it sufficiently converges. 


D. A POSTERIORI VARIANCE OF UNIT WEIGHT 
The sum of the squares of the residuals V7 P V are found by, 


VTPV=LTPL+XTATPL 
The a posteriori variance of unit weight may be calculated, 


VTPV 
n— u 


A2 
Co 


where, n — u is the degree of freedom. 


E. VARIANCE-COVARIANCE MATRIX OF ADJUSTED PARAMETERS 


The variance-covariance matrix for the adjusted parameters is given by, 
Es =(A 7 PAN AT PELPA(A PAY” 
but with P = 032;! this reduces to, 
Es=0 (A PA)” 
F. VARIANCE-COVARIANCE MATRIX OF THE ADJUSTED, OBSERVED 


QUANTITIES 
This matrix is given by, 


Er = ag A(A PA) A 


G. SUMMARY 
To summarize the theoretical aspects of the least squares method using the nonlin- 


ear, observation equation model, 


e The nonlinear system of equations is linearized by expanding in a Taylor series and 
nonlinear terms are truncated (Note: the number of equations must be greater than 
or equal to the number of parameters), 
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e A suitable noise model is chosen as represented by the observation variance- 


covariance matrix and a priori variance of unit weight, and 


e the sum of the squares of the residuals are minimized to yield a minimum variance 


solution 


The operational procedures to estimate values for the unknown parameters are, 


Develop a mathematical mode] that expresses the observations as a function of 
some unknown parameters, 


select appropriate approximate values for the parameters X° (e.g., the receiver co- 
ordinates), 


Form the vector of computed observations L° by evaluating the mathematical 
model at M, 


Determine L, which is the difference vector of the computed observations minus 
the actual observations, 


Form the weight matrix P from the observation variance-covariance matrix and the 
a priori variance of unit weight 


take the partial derivatives of the mathematical model with respect to the un- 
knowns and evalute at X° to form the elements of the design matrix A, 


solve for the change in parameters, X, 


Compute the parameter variance-covariance matrix, a posteriori variance, the resi- 
dual vector, etc., 


apply X to the initial approximations and stop if X has sufficiently converged, 
otherwise repeat the process using X* as the new X° (i.e., iterate) 
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APPENDIX C. CONSTANTS USED IN GASP 


Table 21 presents the constants used in the GASP program. The constants are 


identified; and the common symbols, the units (if applicable), and the values are given. 


Table 15. CONSTANTS USED IN GASP 


Symbol 
WGS-84 Earth’s Semi-major Axis A ee 
WGS-84 Inverse Earth Flattening EMS es 


[Earth's Rotation Raw | o | rae 
A 
[TI 4100 Receiver Li Frequency Ofer | | m 
[T4100 Receiver L2 Frequency offer | | m | 


100 


Value 
6378.137 
298.257223563 


3.986005 x 105 


7.2921151467 x 10- 
3.1415926535897932 
299,792.458 
10.23 x 105 
154 f 
120 fo 
19.029 
24.421 
— 6000 
7600 


APPENDIX D. TABLES OF RESULTS IN CT COORDINATES 


A. BEACH LAB RESULTS 

The results previously expressed in terms of local east, north, and up coordinates 
are presented here as component differences in conventional terrestrial coordinates (X, 
Y, and Z) for the reader more comfortable with this reference frame. Tables 15 through 


17 display these results. No corresponding target plots are shown. 


Table 16. BEACH LAB COMPONENT DIFFERENCES (AX, AY, AZ) USING BROADCAST 
EPHEMERIS: Solution Differences Between GPS and Transit Doppler Methods, 
Second Term is GASP Estimated Formal Error 


14100 
338 










341 
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Table 17. BEACH LAB COMPONENT DIFFERENCES (AX, AY, AZ) USING PRECI 
EPHEMERIS: Solution Differences Between GPS and Transit Doppler Methc 
Second Term is GASP Estimated Formal Error 


DFE 
338 ~09 
“024 4 
“044 4 
339 “0544 
“094 
“1.04 4 
340 EY 
“064 
E 
sa 0554 
0554 
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Table 18. BEACH LAB REPEATABILITY (AX, AY, AZ): AVERAGE OVER DAYS BY RE- 
CEIVER 


Mean AX (m) Mean AY (m) | Mean AZ (m) 
i ee... 


TI A100 “06412 
Broadcast 15 24 IES 
LTES G 


AAA AAA... AAA AAA | 
TI 4100 0.7 + 1.3 0.4 + 1.6 0.1 + 1.1 
Precise 








07412 
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B. LOBOS3 RESULTS 
Tables 18 and 19 display the LOBOS3 results in conventional terrestrial X, Y, and 


Z coordinates. Associated target plots are not presented. 


Table 19. LOBOS3 COMPONENT DIFFERENCES (AX, AY, AZ): 


Solution Differences 


tween GPS and Transit Doppler Methods for Broadcast and Precise Ephemerides, § 
ond Term is GASP Estimated Formal Error 


Day of Year 1991 AX (m) AY (m) 


Broadcast 


Precise 


“06+ 5 “15+ 6 


OE 
1057 
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AZ (m) 


—0.2 + .4 
—0.8 + .5 
0.8 + .4 


0.0 + .4 
—0.5 +.3 
0.2 + .4 


Table 20. LOBOS3 REPEATABILITY (AX, AY, AZ): AVERAGE 
OVER DAYS 


Mean AX (m) | Mean AY(m) | Mean AZ (m) 


D | 
—1.0 + 12 —5.5 + 2.5 
0.4 +.5 









105 


C. BLDG 224.3 RESULTS 
Tab.. 20 presents the solution differences for BLDG 224.3 in the CT coordinate 
system. 


Table 21. BLDG 224.3 COMPONENT DIFFERENCES (AX, AY, AZ): Selective Availab 
Test - GPS Solution Differences Between Previously Determined Coordinates and 
ordinates Determined After Activation of S/A (Using Broadcast Ephemeris) 


CA |  -—297+13 —4.5 + 1.1 5.8 + .6 
ART estes | ero | nn 
CN ee VES 
© D | -—114+10 —8.1 + 1.3 12.4 + .7 
—7.5 + 1.0 2.9 + 1.1 0.7 +.5 
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APPENDIX E. SUMMARIES OF POSITIONING RESULTS 


A. SUMMARY OF TRANSIT POSITION RESULTS 


Station name: DOPS 
Method: Transit Doppler Observations (@DOP3) 
and Conventional Terrestrial Surveying 


Period of Occupation: January 26-30, 1991 

Transit Receiver: Magnavox MX-1502 

Satellites Observed: 77,105,115,124,126,128 

Passes Accepted: 57 

RMS of Residuals: 8.0 cm 

WGS 84 Estimated Station Coordinates 

Monument (X,Y,Z) (meters): -2707273.74 -4353292.80 3781989.36 
Geodetic Coordinates (monument): 


® = 36°36'06.2347 A = 238°07'22.404 h= —23.55 meters 


Standard Deviations (D,A,h) (meters): 
1.0 1.0 1.0 -- 1 Sigma 
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Station name: DOPI 

Method: Transit Doppler Observations (@DOP3) 
and Conventional Terrestrial Surveying 

Period of Occupation: January 26-30, 1991 

Transit Receiver: Magnavox MX-1502 

Satellites Observed: 77,105,115,124,126,128 

Passes Accepted: Sy 

RMS of Residuals: 8.0 cm 

WGS 84 Estimated Station Coordinates 

Monument (X,Y,Z) (meters): -2707255.23 -4353301.05 3781993.34 

Geodetic Coordinates (monument): 


Dd = 36°36'06.3917 A= 238°07'23.211 h = —23.40 meters 


Standard Deviations (,A,h) (meters): 
10 1.0 1.0 -- 1 Sigma 


Station name: DOP2 

Method: Transit Doppler Observations (@ DOP3) 
and Conventional Terrestrial Surveying 

Period of Occupation: January 26-30, 1991 

Transit Receiver: Magnavox MX-1502 

Satellites Observed: 77,105,115,124,126,128 

Passes Accepted: 57 

RMS of Residuals: 8.0 cm 

WGS 84 Estimated Station Coordinates 

Monument (X,Y,Z) (meters): -2707259.43 -4353297.54 3781994.29 

Geodetic Coordinates (monument): 


® = 36°36'06.4311 A =238°07'22.9936 h=—23.45 meters 


Standard Deviations (D,A,h) (meters): 
1.0 1.0 1.0 -- 1 Sigma 
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Station name: LOBOS3 

Method: Transit Doppler Observations 
Period of Occupation: : January 31- February 4, 1991 
Transit Receiver: Magnavox MX-1502 
Satellites Observed: 77,105,115,124,126,128 

Passes Accepted: 62 

RMS of Residuals: 9.0 cm 


WGS 84 Estimated Station Coordinates 


Monument (X,Y,Z) (meters): -2693399.38 -4346211.98 3799864.81 


Geodetic Coordinates (monument): 


D = 36°48'09.683 A = 238°12'46.841 h= —30.45 meters 


Standard Deviations (®,A,h) (meters): 
1.0 1.0 1.0 -- 1 Sigma 
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B. SUMMARY OF GPS POSITION RESULTS 
1. Beach Lab Sites 


Station name: DOPS 
Method: GPS WGS 84 Absolute Point Positioning 
Software: STARPREP (version 1.1), GASP (version 2.0 
Ephemerides: Broadcast 
Satellite Clock States: Broadcast 
Date of Occupation: Day 338, 1990 
GPS Receiver: TI 4100 BEPP/CORE (versions 3.7/5.1) 
Date of Point Position Estimation: June 1991 
Data Collection Span: 6.2 Hours 
PRN Numbers Tracked: 2,6,9,11,12,13,18,19 
Final Number of GASP Observables: 1030 
Percentage of Data Rejected: 2.2 % 
RMS of Residuals: 4.709 cm 
WGS 84 Estimated Station Coordinates 
Monument (X,Y,Z) (meters): -2707274.8418 -4353293.3659 3781988.4006 
Standard Deviations (X,Y,Z) (meters): 
0.6310 0.5967 0.4054 


Geodetic Coordinates (monument): 


D = 36°36'06.1889 A = 238°07'22.3782 h=-—23.268 meters 
Antenna height (monument to electrical center of antenna): 


1.5403 meters 
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Station name: DOPI 
Method: GPS WGS 84 Absolute Point Positioning 
Software: STARPREP (version 1.1), GASP (version 2.0 
Ephemerides: Broadcast 
Satellite Clock States: Broadcast 
Date of Occupation: Day 339, 1990 
GPS Receiver: TI 4100 BEPP/CORE (versions 3.7/5.1) 
Date of Point Position Estimation: June 1991 
Data Collection Span: 6.2 Hours 
PRN Numbers Tracked: Z,6,9,11,12,13, 18,19 
Final Number of GASP Observables: 1017 
Percentage of Data Rejected: 1.2 % 
RMS of Residuals: 4.794 cm 
WGS 84 Estimated Station Coordinates 
Monument (X,Y,Z) (meters): -2707256.7432 -4353302.7647 3781993.1333 
Standard Deviations (X,Y,Z) (meters): 
0.6480 0.6657 0.4016 
Geodetic Coordinates (monument): 


Dd = 36°36'06.3426 A = 238°07'23.1963 A= —21.712 meters 


Antenna height (monument to electrical center of antenna): 


1.3673 meters 
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Station name: DOPS 
Method: GPS WGS 84 Absolute Point Positioning 
Software: STARPREP (version 1.1), GASP (version 2.0 
Ephemerides: Broadcast 
Satellite Clock States: Broadcast 
Date of Occupation: Day 340, 1990 
GPS Receiver: TI 4100 BEPP/CORE (versions 3.7/5.1) 
Date of Point Position Estimation: June 1991 
Data Collection Span: 7.1 Hours 
PRN Numbers Tracked: 2,6,9,11,12,13,18,19 
Final Number of GASP Observables: 1185 
Percentage of Data Rejected: 2.3 % 
RMS of Residuals: 4.584 cm 
WGS 84 Estimated Station Coordinates 
Monument (X,Y,Z) (meters): -2707276.6237 -4353293.1626 3781989.1370 
Standard Deviations (X,Y,Z) (meters): 
0.5496 0.5643 0.3848 
Geodetic Coordinates (monument): 


® = 36°36'06.1932 A = 238°07'22.3130 A= —22.212 meters 


Antenna height (monument to electrical center of antenna): 


1.5873 meters 





112 


Station name: DOPS 
Method: GPS WGS 84 Absolute Point Positioning 
Software: STARPREP (version 1.1), GASP (version 2.0 
Ephemerides: Broadcast 
Satellite Clock States: Broadcast 
Date of Occupation: Day 341, 1990 
GPS Receiver: TI 4100 BEPP/CORE (versions 3.7/5.1) 
Date of Point Position Estimation: June 1991 
Data Collection Span: 6.5 Hours 
PRN Numbers Tracked: 2,6,9,11,12,13,18,19 
Final Number of GASP Observables: 690 
Percentage of Data Rejected: 38.2 % 
RMS of Residuals: 3.306 cm 
WGS 84 Estimated Station Coordinates 
Monument (X,Y,Z) (meters): -2107272.5806 -4353293.6093 3781988.4578 
Standard Deviations (X,Y,Z) (meters): 
0.5463 0.5307 0.3803 
Geodetic Coordinates (monument): 


Dd = 36°36'06.2095 A = 238°07'22.4606 h = —24.027 meters 


Antenna height (monument to electrical center of antenna): 


1.5103 meters 
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Station name: DOPI 
Method: GPS WGS 84 Absolute Point Positioning 
Software: STARPREP (version 1.1), GASP (version 2.0 
Ephemerides: Broadcast 
Satellite Clock States: Broadcast 
Date of Occupation: Day 338, 1990 
GPS Receiver: Ashtech LD XII 
Receiver software version 
Date of Point Position Estimation: June 1991 
Data Collection Span: 5.6 Hours 
PRN Numbers Tracked: 2,6,9,11,12,15,16,18,19 
Final Number of GASP Observables: 1733 
Percentage of Data Rejected: 9,3 % 
RMS of Residuals: 3.336 cm 
WGS 84 Estimated Station Coordinates 
Monument (X,Y,Z) (meters): -2707254.1323 -4353299.9074 3781992.3652 
Standard Deviations (X,Y,Z) (meters): 
0.6568 0.4494 0.4107 


Geodetic Coordinates (monument): 


D = 36°36'06.3962 A = 238°07'23.2248 h = —25.225 meters 
Antenna height (monument to electrical center of antenna): 


1.3313 meters 
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Station name: DOP5 
Method: GPS WGS 84 Absolute Point Positioning 
Software: STARPREP (version 1.1), GASP (version 2.0 
Ephemendes: Broadcast 
Satellite Clock States: Broadcast 
Date of Occupation: Day 339, 1990 
GPS Receiver: Ashtech LD XII 
Receiver software version 

Date of Point Position Estimation: June 1991 
Data Collection Span: 6 Hours 
PRN Numbers Tracked: 2,6,9,11,12,15,16,18,19 
Final Number of GASP Observables: 1079 
Percentage of Data Rejected: 19.8% 
RMS of Residuals: 4.590 cm 
WGS 84 Estimated Station Coordinates 
Monument (X,Y,Z) (meters): -2707269.1762 -4353291.6786 3781988.0816 
Standard Deviations (X,Y,Z) (meters): 

| 0.7447 0.4817 0.4241 


Geodetic Coordinates (monument): 


Dd = 36°36'06.2662 A =238°07'22.5359 h=-—27.011 meters 


Antenna height (monument to electrical center of antenna): 


1.3373 meters 
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Station name: DOPI 
Method: GPS WGS 84 Absolute Point Positioning 
Software: STARPREP (version 1.1), GASP (version 2.0 
Ephemerides: Broadcast 
Satellite Clock States: Broadcast 
Date of Occupation: Day 340, 1990 
GPS Receiver: Ashtech LD XII 
Receiver software version 
Date of Point Position Estimation: June 1991 
Data Collection Span: 6.8 Hours 
PRN Numbers Tracked: 2,6,9,11,12,15,16,18,19 
Final Number of GASP Observables: 1297 
Percentage of Data Rejected: 9.7 % 
RMS of Residuals: 4.481 cm 
WGS 84 Estimated Station Coordinates 
Monument (X,Y,Z) (meters): -2707253.7859 -4353300.2267 3781992.1104 
Standard Deviations (X,Y,Z) (meters): 
0.6444 0.4487 0.4048 


Geodetic Coordinates (monument): 


Dd = 36°36'06.3879 A= 238°07'23.2434 h=—25.306 meters 
Antenna height (monument to electrical center of antenna): 


1.2483 meters 





116 


Station name: DOP! 
Method: GPS WGS 84 Absolute Point Positioning 
Software: STARPREP (version 1.1), GASP (version 2.0 
Ephemerides: Broadcast 
Satellite Clock States: Broadcast 
Date of Occupation: Day 341, 1990 
GPS Receiver: Ashtech LD XII 
Receiver software version 
Date of Point Position Estimation: June 1991 
Data Collection Span: 6.9 Hours 
PRN Numbers Tracked: 2,6,9,11,12,15,16,18,19 
Final Number of GASP Observables: 1255 


Percentage of Data Rejected: 13.6% 
RMS of Residuals: 4.259 cm 
WGS 84 Estimated Station Coordinates 
Monument (X,Y,Z) (meters): -2707256.4461 -4353299.5540 3781993.3044 
Standard Deviations (X,Y,Z) (meters): 
0.6027 0.4345 0.3941 


Geodetic Coordinates (monument): 
® = 36°36'06.4029 A = 238°07'23.1382 h=—23.925 meters 
Antenna height (monument to electrical center of antenna): 


1.4823 meters 
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Station name: DOP2 
Method: GPS WGS 84 Absolute Point Positioning 
Software: STARPREP (version 1.1), GASP (version 2.0 
Ephemerides: Broadcast 
Satellite Clock States: Broadcast 
Date of Occupation: Day 338, 1990 
GPS Receiver: Trimble 4000ST 
Receiver software version 4.3X 
Date of Point Position Estimation: June 1991 
Data Collection Span: 6 Hours 
PRN Numbers Tracked: 2,6,9,11,12,15,16,18,19 
Final Number of GASP Observables: 1151 
Percentage of Data Rejected: 11.1% 
RMS of Residuals: 4.014 cm 
WGS 84 Estimated Station Coordinates 
Monument (X,Y,Z) (meters): -2707258.5156 -4353294.9312 3781993.6745 
Standard Deviations (X,Y,Z) (meters): 
0.6452 0.4227 0.3927 
Geodetic Coordinates (monument): 


D = 36°36'06.4673 A = 238°07'22.9693 h = —25.978 meters 


Antenna height (monument to electrical center of antenna): 


1.3313 meters 
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Station name: DOP2 
Method: GPS WGS 84 Absolute Point Positioning 
Software: STARPREP (version 1.1), GASP (version 2.0 
Ephemerides: Broadcast 
Satellite Clock States: Broadcast 
Date of Occupation: Day 339, 1990 
GPS Receiver: Trimble 4000ST 
Receiver software version 4.3X 
Date of Point Position Estimation: June 1991 
Data Collection Span: 6 Hours 
PRN Numbers Tracked: 2,6,9,11,12,15,16,18,19 
Final Number of GASP Observables: 1073 
Percentage of Data Rejected: 18.0% 
RMS of Residuals: 3.929 cm 
WGS 84 Estimated Station Coordinates 
Monument (X,Y,Z) (meters): -2707254.9181 -4353295.7659 3781992.9337 
Standard Deviations (X,Y,Z) (meters): 
0.6625 0.4241 0.3944 


Geodetic Coordinates (monument): 


® = 36°36'06.4710 A = 238°07'23.1099 h= —27.376 meters 
Antenna height (monument to electrical center of antenna): 


1.3023 meters 
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Station name: 
Method: 

Software: 
Ephemerides: 
Satellite Clock States: 
Date of Occupation: 
GPS Receiver: 


Date of Point Position Estimation: 
Data Collection Span: 

PRN Numbers Tracked: 

Final Number of GASP Observables: 


Percentage of Data Rejected: 
RMS of Residuals: 
WGS 84 Estimated Station Coordinates 


Monument (X,Y,Z) (meters): 


Standard Deviations (X,Y,Z) (meters): 


Geodetic Coordinates (monument): 


DOP2 

GPS WGS 84 Absolute Point Positioning 
STARPREP (version 1.1), GASP (version 2.0 
Broadcast 

Broadcast 

Day 340, 1990 

Trimble 4000ST 

Receiver software version 4.3X 

June 1991 

7.5 Hours 

2,6,9,11,12,15,16,18,19 

1282 

18.8% 

3.896 cm 


-2707258.2384 -4353296.3532 3781994.0224 


0.5624 0.4229 0.3724 


D = 36°36'06.4558 A = 238°07'23.0090 h=—24.919 meters 


Antenna height (monument to electrical center of antenna): 


1.2743 meters 
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Station name: DOP2 
Method: GPS WGS 84 Absolute Point Positioning 
Software: STARPREP (version 1.1), GASP (version 2.0 
Ephemerides: Broadcast 
Satellite Clock States: Broadcast 
Date of Occupation: Day 341, 1990 
GPS Receiver: Trimble 4000ST 
Receiver software version 4.3X 
Date of Point Position Estimation: June 1991 
Data Collection Span: 6.7 Hours 
PRN Numbers Tracked: 2,6,9,11,12,15,16,18,19 
Final Number of GASP Observables: 1280 


Percentage of Data Rejected: 13.5% 
RMS of Residuals: 4.096 cm 
WGS 84 Estimated Station Coordinates 
Monument (X,Y,Z) (meters): -2707259.2464 -4353297.8957 3781994.4146 
Standard Deviations (X,Y,Z) (meters): 
0.6003 0.4309 0.3825 


Geodetic Coordinates (monument): 
D = 36°36'06.4304 A = 238°07'23.0073 h = —23.206 meters 
Antenna height (monument to electrical center of antenna): 


1.3933 meters 
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Station name: DOPS 
Method: GPS WGS 84 Absolute Point Positioning 
Software: STARPREP (version 1.1), GASP (version 2.0 
Ephemerides: Precise WGS 84 EF##90336 
Satellite Clock States: Precise PC90336 
Date of Occupation: Day 338, 1990 
GPS Receiver: TI 4100 BEPP/CORE (versions 3.7/5.1) 
Date of Point Position Estimation: June 1991 
Data Collection Span: 6.2 Hours 
PRN Numbers Tracked: 2,6,9,11,12,13,18,19 
Final Number of GASP Observables: 1033 
Percentage of Data Rejected: 1.9 % 
RMS of Residuals: 4.717 cm 
WGS 84 Estimated Station Coordinates 
Monument (X,Y,Z) (meters): -2707271.5459 -4353291.7286 3781989.5013 
Standard Deviations (X,Y,Z) (meters): 
0.6306 0.5974 0.4060 


Geodetic Coordinates (monument): 


® = 36°36'06.2781 A = 238°07'22.4560 h= —25.125 meters 
Antenna height (monument to electrical center of antenna): 


1.5403 meters 
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Station name: DOPl 

Method: GPS WGS 84 Absolute Point Positioning 
Software: STARPREP (version 1.1), GASP (version 2.0 
Ephemerides: Precise WGS84 EF##90336 

Satellite Clock States: Precise PC90336 

Date of Occupation: Day 339, 1990 

GPS Receiver: TI 4100 BEPP/CORE (versions 3.7/5.1) 
Date of Point Position Estimation: June 1991 

Data Collection Span: 6.2 Hours 

PRN Numbers Tracked: 2,6,9,11,12,13,18,19 

Final Number of GASP Observables: 1021 


Percentage of Data Rejected: 0.8 % 
RMS of Residuals: 4.361 cm 
WGS 84 Estimated Station Coordinates 
Monument (X,Y,Z) (meters): -2707255.4509 -4353300.4770 3781992.9055 
Standard Deviations (X,Y,Z) (meters): 
0.5886 0.6084 0.3707 


Geodetic Coordinates (monument): 
® = 36°36'06.3875 A =238°07'23.1918 h= —23.955 meters 
| Antenna height (monument to electrical center of antenna): 


| iea meters 
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Station name: DOPS 
Method: GPS WGS 84 Absolute Point Positioning 
Software: STARPREP (version 1.1), GASP (version 2.0 
Ephemerides: Precise WGS 84 EF##90336 
Satellite Clock States: Precise PC90336 
Date of Occupation: Day 340, 1990 
GPS Receiver: TI 4100 BEPP/CORE (versions 3.7/5.1) 
Date of Point Position Estimation: June 1991 
Data Collection Span: 7.1 Hours 
PRN Numbers Tracked: 2,6,9,11,12,13,18,19 
Final Number of GASP Observables: 1182 
Percentage of Data Rejected: 2.6 % 
RMS of Residuals: 4.563 cm 
WGS 84 Estimated Station Coordinates 
Monument (X,Y,Z) (meters): -2707274.3957 -4353290.8606 3781988.3570 
Standard Deviations (X,Y,Z) (meters): 
0.5548 0.5665 0.3870 
Geodetic Coordinates (monument): 


D = 36°36'06.2334 A =238°07'22.3402 h = —25.191 meters 


Antenna height (monument to electrical center of antenna): 


1.5873 meters 
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Station name: DOPS 
Method: GPS WGS 84 Absolute Point Positioning 
Software: STARPREP (version 1.1), GASP (version 2.0 
Ephemerides: Precise WGS 84 EF##90336 
Satellite Clock States: Precise PC90336 
Date of Occupation: Day 341, 1990 
GPS Receiver: TI 4100 BEPP/CORE (versions 3.7/5.1) 
Date of Point Position Estimation: June 1991 
Data Collection Span: 6.5 Hours 
PRN Numbers Tracked: 2,6,9,11,12,13,18,19 
Final Number of GASP Observables: 1090 
Percentage of Data Rejected: 2.4 % 
RMS of Residuals: 4.506 cm 
WGS 84 Estimated Station Coordinates 
Monument (X,Y,Z) (meters): -2707272.3968 -4353294.6322 3781990.9281 
Standard Deviations (X,Y,Z) (meters): 
0.5813 0.5608 0.3847 
Geodetic Coordinates (monument): 


D = 36°36'06.2589 A = 238°07'22.4886 A= —21.934 meters 


Antenna height (monument to electrical center of antenna): 


1.5103 meters 
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Station name: DOPI 
Method: GPS WGS 84 Absolute Point Positioning 
Software: STARPREP (version 1.1), GASP (version 2.0 
Ephemerides: Precise WGS 84 EF##90336 
Satellite Clock States: Precise PC90336 
Date of Occupation: Day 338, 1990 
GPS Receiver: Ashtech LD XII 
Receiver software version 
Date of Point Position Estimation: June 1991 
Data Collection Span: 5.6 Hours 
PRN Numbers Tracked: 2,6,9,11,12,15,16,18,19 
Final Number of GASP Observables: 1737 
Percentage of Data Rejected: 9.1% 
RMS of Residuals: 3.265 cm 
WGS 84 Estimated Station Coordinates 
Monument (X,Y,Z) (meters): -2707256.1438 -4353300.1852 3781992.4439 
Standard Deviations (X,Y,Z) (meters): 
0.6395 0.4393 0.4035 


Geodetic Coordinates (monument): 


®D = 36°36'06.3732 A= 238°07'23.1619 h= —24.136 meters 
Antenna height (monument to electrical center of antenna): 


1.3313 meters 
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Station name: DOPS 
Method: GPS WGS 84 Absolute Point Positioning 
Software: STARPREP (version 1.1), GASP (version 2.0 
Ephemerides: Precise WGS 84 EF##90336 
Satellite Clock States: Precies PC90336 
Date of Occupation: Day 339, 1990 
GPS Receiver: Ashtech LD XII 
Receiver software version 
Date of Point Position Estimation: June 199] 
Data Collection Span: 6 Hours 
PRN Numbers Tracked: 2,6,9,11,12,15,16,18,19 
Final Number of GASP Observables: 1092 
Percentage of Data Rejected: 18.8% 
RMS of Residuals: 4.347 cm 
WGS 84 Estimated Station Coordinates 
Monument (X,Y,Z) (meters): -2707274.2183 -4353292.9093 3781988.8495 
Standard Deviations (X,Y,Z) (meters): 
0.7005 0.4521 0.4073 
Geodetic Coordinates (monument): 


® = 36°36'06.2144 A = 238°07'22.3898 h= —23.576 meters 


Antenna height (monument to electrical center of antenna): 


1.3373 meters 
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Station name: DOPI 
Method: GPS WGS 84 Absolute Point Positioning 
Software: STARPREP (version 1.1), GASP (version 2.0 
Ephemerides: Precise WGS 84 EF##90336 
Satellite Clock States: Precise PC90336 
Date of Occupation: Day 340, 1990 
GPS Receiver: Ashtech LD XII 
Receiver software version 
Date of Point Position Estimation: June 1991 
Data Collection Span: 6.8 Hours 
PRN Numbers Tracked: 2,6,9,11,12,15,16,18,19 
Final Number of GASP Observables: 1297 
Percentage of Data Rejected: 9.7% 
RMS of Residuals: 4.338 cm 
WGS 84 Estimated Station Coordinates 
Monument (X,Y,Z) (meters): -2707254.4187 -4353299.6465 3781991.6860 
Standard Deviations (X,Y,Z) (meters): 
0.6276 0.4352 0.3971 
Geodetic Coordinates (monument): 


® = 36°36'06.3799 A =238°07'23.2094 h = —25.686 meters 


Antenna height (monument to electrical center of antenna): 


1.2483 meters 
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Station name: DOP! 
Method: GPS WGS 84 Absolute Point Positioning 
Software: STARPREP (version 1.1), GASP (version 2.0 
Ephemerides: Precise WGS 84 EF##90336 
Satellite Clock States: Precise PC90336 
Date of Occupation: Day 341, 1990 
GPS Receiver: Ashtech LD XI] 
Receiver software version 
Date of Point Position Estimation: June 1991 
Data Collection Span: 6.9 Hours 
PRN Numbers Tracked: 2,6,9,11,12,15,16,18,19 
Final Number of GASP Observables: 1263 


Percentage of Data Rejected: 13.1% 
RMS of Residuals: 4.420 cm 


WGS 84 Estimated Station Coordinates 
Monument (X,Y,Z) (meters): -2707255.1764 -4353300.5816 3781993.8009 
Standard Deviations (X,Y,Z) (meters): 
l 0.6308 0.4529 0.4070 
Geodetic Coordinates (monument): 


® = 36°36'06.4119 A = 238°07'23.2034 h = —23.466 meters 
Antenna height (monument to electrical center of antenna): 


1.4823 meters 





129 


Station name: DOP2 
Method: GPS WGS 84 Absolute Point Positioning 
Software: STARPREP (version 1.1), GASP (version 2.0 
Ephemerides: Precise WGS 84 EF##90336 
Satellite Clock States: Precise PC90336 
Date of Occupation: Day 338, 1990 
GPS Receiver: Trimble 4000ST 
Receiver software version 4.3X 
Date of Point Position Estimation: June 1991 
Data Collection Span: 6 Hours 
PRN Numbers Tracked: 2,6,9,11,12,15,16,18,19 
Final Number of GASP Observables: 1164 
Percentage of Data Rejected: 11.4% 
RMS of Residuals: 3.807 cm 
WGS 84 Estimated Station Coordinates 
Monument (X,Y,Z) (meters): -2707260.5913 -4353295.8412 3781994.0927 
Standard Deviations (X,Y,Z) (meters): 
0.6084 0.3998 0.3763 
Geodetic Coordinates (monument): 


D = 36°36'06.4420 A =238°07'22.9177 h=—24.228 meters 


Antenna height (monument to electrical center of antenna): 


1.3313 meters 
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Station name: DOP2 
Method: GPS WGS 84 Absolute Point Positioning 
Software: STARPREP (version 1.1), GASP (version 2.0 
Ephemerides: Precise WGS 84 EF##90336 
Satellite Clock States: Precise PC90336 
Date of Occupation: Day 339, 1990 
GPS Receiver: Trimble 4000ST 
Receiver software version 4.3X 
Date of Point Position Estimation: June 1991 
Data Collection Span: 6 Hours 
PRN Numbers Tracked: 2,6,9,11,12,15,16,18,19 
Final Number of GASP Observables: 1089 


Percentage of Data Rejected: 16.8% 
RMS of Residuals: 3.715 cm 
WGS 84 Estimated Station Coordinates 
Monument (X,Y,Z) (meters): -2707258.8586 -4353296.6201 3781993.4158 
Standard Deviations (X,Y,Z) (meters): 
0.6186 0.3956 0.3804 


Geodetic Coordinates (monument): 
D = 36°36'06.4293 A= 238°07'22.9934 h=—24.836 meters 
Antenna height (monument to electrical center of antenna): 


1.3023 meters 
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Station name: DOP2 
Method: GPS WGS 84 Absolute Point Positioning 
Software: STARPREP (version 1.1), GASP (version 2.0 
Ephemerides: Precise WGS 84 EF##90336 
Satellite Clock States: Precise PC90336 
Date of Occupation: Day 340, 1990 
GPS Receiver: Trimble 4000ST 
Receiver software version 4.3X 
Date of Point Position Estimation: June 1991 
Data Collection Span: 7.5 Hours 
PRN Numbers Tracked: 2,6,9,11,12,15,16,18,19 
Final Number of GASP Observables: 1281 
Percentage of Data Rejected: 18.8% 
RMS of Residuals: 3.793 cm 
WGS 84 Estimated Station Coordinates 
Monument (X,Y,Z) (meters): -2707259.1230 -4353296.3184 3781993.6501 
Standard Deviations (X,Y,Z) (meters): 
0.5582 0.4141 0.3703 


Geodetic Coordinates (monument): 


Dd = 36°36'06.4377 A= 238°07'22.9780 h = —24.789 meters 
Antenna height (monument to electrical center of antenna): 


1.2743 meters 
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Station name: DOP2 
Method: GPS WGS 84 Absolute Point Positioning 
Software: STARPREP (version 1.1), GASP (version 2.0 
Ephemerides: Precise WGS 84 EF##90336 
Satellite Clock States: Precise PC90336 
Date of Occupation: Day 341, 1990 
GPS Receiver: Trimble 4000ST 
Receiver software version 4.3X 
Date of Point Position Estimation: June 1991 
Data Collection Span: 6.7 Hours 
PRN Numbers Tracked: 2,6,9,11,12,15,16,18,19 
Final Number of GASP Observables: 1285 
Percentage of Data Rejected: 13.1% 
RMS of Residuals: 4.214 cm 
WGS 84 Estimated Station Coordinates 
Monument (X,Y,Z) (meters): -2707258.4303 -4353298.5061 3781994.7611 
Standard Deviations (X,Y,Z) (meters): 
0.6086 0.4354 0.3854 
Geodetic Coordinates (monument): 


® = 36°36'06.4377 A = 238°07'23.0481 h = —22.929 meters 


Antenna height (monument to electrical center of antenna): 


1.3933 meters 
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2. LOBOS3 Site 


Station name: LOBOS3 
Method: GPS WGS 84 Absolute Point Positioning 
Software: STARPREP (version 1.1), GASP (version 2.0 
Ephemerides: Broadcast 
Satellite Clock States: Broadcast 
Date of Occupation: Day 36, 1991 
GPS Receiver: TI 4100 BEPP/CORE (versions 3.7/5.1) 
Date of Point Position Estimation: June 1991 
Data Collection Span: 8.6 Hours 
PRN Numbers Tracked: 3,6,11,12,13,16,17,20 
Final Number of GASP Observables: 964 
Percentage of Data Rejected: 2.0 % 
RMS of Residuals: 4.847 cm 
WGS 84 Estimated Station Coordinates 
Monument (X,Y,Z) (meters): -2693399,4335 -4346214.6434 3799864.5746 
Standard Deviations (X,Y,Z) (meters): 
0.6637 0.5214 0.3612 
Geodetic Coordinates (monument): 


Dd = 36°48’09.6322 A =238°12'46.8955 h = —28.749 meters 


Antenna height (monument to electrical center of antenna): 


1.5673 meters 
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Station name: LOBOS3 
Method: GPS WGS 84 Absolute Point Positioning 
Software: STARPREP (version 1.1), GASP (version 2.0 
Ephemerides: Broadcast 
Satellite Clock States: Broadcast 
Date of Occupation: Day 37, 1991 
GPS Receiver: TI 4100 BEPP/CORE (versions 3.7/5.1 
Date of Point Position Estimation: June 1991 i 
Data Collection Span: 10.8 Hours 
PRN Numbers Tracked: 3,6,11,12,13,16,17,20 
Final Number of GASP Observables: 1361 
Percentage of Data Rejected: 0.9 % 
RMS of Residuals: 6.786 cm 
WGS 84 Estimated Station Coordinates 
Monument (X,Y,Z) (meters): -2693400.0239 -4346219.4479 3799864.0505 
Standard Deviations (X,Y,Z) (meters): 

0.8974 0.6011 0.4497 


Geodetic Coordinates (monument): 


Hd = 36°48'09.5332 A =238°12'46.9773 h=—25.544 meters 


Antenna height (monument to electrical center of antenna): 


1.5673 meters 
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Station name: LOBOS3 
Method: GPS WGS 84 Absolute Point Positioning 
Software: STARPREP (version 1.1), GASP (version 2.0 
Ephemerides: Broadcast 
Satellite Clock States: Broadcast 
Date of Occupation: Day 38, 1991 
GPS Receiver: TI 4100 BEPP/CORE (versions 3.7/5.1) 
Date of Point Position Estimation: June 1991 
Data Collection Span: 10.8 Hours 
PRN Numbers Tracked: 3,6,11,12,13,16,17,20 
Final Number of GASP Observables: 1293 
Percentage of Data Rejected: 0.8 % 
RMS of Residuals: 5.733 cm 
WGS 84 Estimated Station Coordinates 
Monument (X,Y,Z) (meters): -2693401.7366 -4346218.3372 3799865.5840 
Standard Deviations (X,Y,Z) (meters): 
0.7575 0.5162 0.3880 


Geodetic Coordinates (monument): 


Dd = 36°48'09.5738 A = 238°12'46.8950 h = —24.659 meters 
Antenna height (monument to electrical center of antenna): 


1.5673 meters 
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Station name: LOBOS3 
Method: GPS WGS 84 Absolute Point Positioning 
Software: STARPREP (version 1.1), GASP (version 2.0 
Ephemerides: Precise WGS 84 EF##91034 
Satellite Clock States: Precise PC91034 
Date of Occupation: Day 36, 1991 
GPS Receiver: TI 4100 BEPP/CORE (versions 3.7/5.1) 
Date of Point Position Estimation: June 1991 
Data Collection Span: 8.6 Hours 
PRN Numbers Tracked: 3,6,11,12,13,16,17,20 
Final Number of GASP Observables: 970 
Percentage of Data Rejected: 1.4% 
RMS of Residuals: 4.947 cm 
WGS 84 Estimated Station Coordinates 
Monument (X,Y,Z) (meters): -2693400.0858 -4346212.4518 3799864.7808 
Standard Deviations (X,Y,Z) (meters): 
0.6707 0.5258 0.3626 
Geodetic Coordinates (monument): 


D = 36°48'09.6671 A = 238°12'46.8266 h = —29.842 meters 


Antenna height (monument to electrical center of antenna): 


1.5673 meters 
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Station name: 

Method: 

Software: 

Ephemerides: 

Satellite Clock States: 

Date of Occupation: 

GPS Receiver: 

Date of Point Position Estimation: 
Data Collection Span: 

PRN Numbers Tracked: 

Final Number of GASP Observables: 
Percentage of Data Rejected: 
RMS of Residuals: 


LOBOS3 

GPS WGS 84 Absolute Point Positioning 
STARPREP (version 1.1), GASP (version 2.0 
Precise WGS 84 EF##91034 

Precise PC91034 

Day 37, 1991 

TI 4100 BEPP/CORE (versions 3.7/5.1) 
June 1991 

10.8 Hours 

3,6,11,12,13,16,17,20 

1364 

0.7% 

4.956 cm 


WGS 84 Estimated Station Coordinates 


Monument (X,Y,Z) (meters): 


Standard Deviations (X,Y,Z) (meters): 


Geodetic Coordinates (monument): 


-2693398.8698 -4346211.9262 3799864.2645 


0.6641 0.4467 0.3459 


D = 36°48'09.6748 A = 238°12'46.8571 h = —31.022 meters 


Antenna height (monument to electrical center of antenna): 


1.5673 meters 
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Station name: LOBOS3 

Method: GPS WGS 84 Absolute Point Positioning 
Software: STARPREP (version 1.1), GASP (version 2.0 
Ephemerides: Precise WGS 84 EF##91034 

Satellite Clock States: Precise PC91034 

Date of Occupation: Day 38, 1991 

GPS Receiver: TI 4100 BEPP/CORE (versions 3.7/5.1) 
Date of Point Position Estimation: June 1991 

Data Collection Span: 10.8 Hours 

PRN Numbers Tracked: 3,6,11,12,13,16,17,20 

Final Number of GASP Observables: 1298 


Percentage of Data Rejected: 0.4% 
RMS of Residuals: 5.160 cm 
WGS 84 Estimated Station Coordinates 
Monument (X,Y,Z) (meters): -2693398.3618 -4346212.7487 3799864.9918 
Standard Deviations (X,Y,Z) (meters): 
0.6893 0.4706 0.3577 


Geodetic Coordinates (monument): 
D = 36°48'09.6858 A = 238°12'46.8920 h= —30.241 meters 
Antenna height (monument to electrical center of antenna): 


1.5673 meters 
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